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I Abstract. Complete, conformally flat metrics of constant positive scalar curvature on 

CO ' the complement of k points in the n-sphere, fc > 2, n > 3, were constructed by R. Schoen 

[S2]. We consider the problem of determining the moduli space of all such metrics. All 
such metrics are asymptotically periodic, and we develop the linear analysis necessary 
to understand the nonlinear problem. This includes a Fredholm theory and asymptotic 
■ regularity theory for the Laplacian on asymptotically periodic manifolds, which is of 

' independent interest. The main result is that the moduli space is a locally real analytic 

variety of dimension fc. For a generic set of nearby conformal classes the moduli space 
is shown to be a fc— dimensional real analytic manifold. The structure as a real analytic 
variety is obtained by writing the space as an intersection of a Fredholm pair of infinite 
dimensional real analytic manifolds. 
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I. Introduction 

Much has been clarified in the past ten years about the behavior of solutions of the 
semilinear elliptic equation relating the scalar curvature functions of two conformally 
related metrics. The starting point for these recent developments was R. Schoen's reso- 
lution of the Yamabe problem on compact manifolds [SI], capping the work of a number 
of mathematicians over many years. Soon thereafter Schoen [S2] and Schoen- Yau [SY] 
made further strides in understanding weak solutions of this equation, particularly on 
the sphere, and its relationship with conformal geometry. Of particular interest here 
is the former, [S2]; in that paper, Schoen constructs metrics with constant positive 
scalar curvature on S'^, conformal to the standard round metric, and with prescribed 
isolated singularities (he also constructs solutions with certain, more general, singular 
sets). This 'singular Yamabe problem' is to find a metric g = u^/^'^~'^^ qq on a domain 
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§'^\A which is complete and has constant scalar curvature R{g). This is equivalent to 
finding a positive function u satisfying 




n-2 



R{g)u 



n + 2 

n-2 — 



on §"\A 



(1.1) 



4(n- 1) 



g complete on §"\A, 



R{g) = constant, 



where R{go) — n(n— 1) is the scalar curvature of the round metric qq. The completeness 
of g requires that u tend to infinity, in an averaged sense, on approach to A. 

The earliest work on this singular Yamabe problem seems to have been that of 
C. Loewner and L. Nirenberg [LN], where metrics with constant negative scalar curva- 
ture are constructed. Later work on this 'negative' case was done by Aviles-McOwen 
[AMc] cf. also [Mc] and Finn-McOwen [FMc] , where the background manifold and met- 
ric are allowed to be arbitrary. For a solution with R{g) < to exist, it is necessary and 
sufiicient that dim(A) > (n — 2)/2. A partial converse is that if a solution with R{g) > 
exists (at least when M — S"^), it is necessary that k < (n — 2)/2, [SY]. Schoen gave the 
first general construction of solutions with R{g) > [S2] . These solutions have singular 
set A which is either discrete or nonrectifiable. Many new solutions on the sphere with 
singular set A, a smooth perturbation of an equatorial fc-sphere, 1 < k < (n — 2)/2, 
are constructed in [MS] . F. Pacard [Pa] has recently constructed positive complete so- 
lutions on the sphere with singular set an arbitrary smooth submanifold of dimension 
(n — 2)/2. Regularity of solutions near the singular set A, no matter the sign of the 
curvature or background manifold or metric, is examined in [Ml], cf. also [ACF] for a 
special case of relevance to general relativity. A more detailed account of part of this 
history is given in [MS] . 

In this paper we return to this problem in the setting studied by Schoen [S2] on 
where A is a finite point-set: 



Hereafter, A will always be taken to be this set, unless indicated specifically, and the 
scalar curvature R{g) attained by the conformal metric will always be n(n — 1). In this 
case, no solutions of (1.1) exist when A; = 1. A proof of this, following from a general 
symmetry theorem, is indicated below. Thus we always assume that k > 2. Since this 
problem is conformally invariant, the set A may be replaced by -F(A) for any conformal 
transformation F G 0{n + 1, 1). A simple topological argument (given in [S2]) shows 
that using such a transformation we can always arrange that A is 'balanced,' i.e. that 
the points pj sum to zero as vectors in R"^"'"-^. Henceforth this will be tacitly assumed as 
well. Notice that if A is balanced, k > 2, and if F is a conformal transformation which 
preserves A, then (since some power of F fixes three points on the sphere, and also the 
origin in R"^+^) F must be an orthogonal transformation. If A; = 2 then F could also 
be a dilation. Also, any balanced configuration A is contained in a minimal equatorial 
subsphere S'^ C and an Alexandrov reflection argument similar to the one in §8 of 
[CGS] shows that if u is an arbitrary solution to (1.1), and F is a rotation preserving 



(1.2) 
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A pointwise, then u must also be invariant under F. Observe that this imphes that no 
complete solutions exist on If u were such a solution, then because any other 

point can be moved to be antipodal to p, this reflection argument would show that u is 
rotationally symmetric with respect to any other point on the sphere. Hence u would 
have to be constant, which contradicts completeness. Hence (except when A; = 2) we do 
not get 'trivial' families of solutions of (1.1) obtained by pulling back a fixed solution 
by a family of conformal transformations. Note, however, that if F preserves A, but 
permutes the points, there is no reason to expect that most solutions of (1.1) will be 
fixed by F; indeed, it is not even clear a priori that F-invariant solutions exist in this 
case. 

Our aim in this work is to consider the moduli space M.a, which is by definition the 
set of all smooth positive functions u to the problem 



where A = {pi, ■ ■ ■ ,pk} is any fixed, balanced set of k points in S"^. We remark that 
the geometric condition that the metric g be complete is equivalent to the analytic one 
of requiring that the set A consist of nonremovable singularities of u. This space might 
be called the 'PDE moduli space' M^^^, to distinguish it from the 'geometric moduli 
space' M.j^°'^ which consists of all geometrically distinct solutions of this problem. 
In most instances these spaces coincide, although the second could conceivably be 
somewhat smaller than the first if some solutions admit nontrivial isometries (e.g. as 
in the case k — 2 discussed in §2 below). A4\ will always denote the former of these 
spaces here. By Schoen's work [S2], Ma is nonempty whenever A; > 1. In fact, his 
construction yields families of solutions. As this was not the aim of his work, this is 
not made explicit there, nor are the free parameters in his construction counted. 

We shall examine a number of questions, both local and global, concerning the nature 
of this moduli space. The simplest of these is whether Ma is a manifold, or otherwise 
tractable set, and if so, what is its dimension. Our main result is 

Theorem 1.4. Ma is locally a real analytic variety of formal dimension k. 

The formal dimension is the dimension predicted by an index theorem. As is usual in 
moduli space theories, obstructions may well exist to prevent Ma from attaining this 
dimension. We clarify this and give a more careful statement in Theorem 5.4, Corollary 
5.5 and Theorem 6.13 below. We also describe how natural parameters on Ma may 
arise. On the linear level, these are given by certain scattering theoretic information 
for the metrics g E Ma- Another geometric description is given by the Pohozaev 
invariants which are defined in §3. We also obtain information on a geometrically 
natural compactification of Ma which is obtained by taking the union with lower- 
dimensional moduli spaces corresponding to singular sets A' C A. We are, as yet, 
unable to provide a satisfactory description of the interior singularities of Al a, or to 



(1.3) 




g = tt"-2 ^0 complete on S'^\A, 



on §"\A 
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determine whether this compactification, AAki is itself a compact real (semi-) analytic 
variety. This latter property is, by all indications, true. We hope to return to this later. 

It is important to mention that there are many similarities between the theory of 
constant scalar curvature metrics on S"\A where A is finite and that of embedded, 
complete, constant mean curvature (CMC) surfaces with k ends in M?. The first exam- 
ples of such CMC surfaces were given by N. Kapouleas [Kap] . One-parameter families 
of solutions with symmetry were constructed by K. Grosse-Brauckmann [B]. Further 
general results on the structure of these surfaces and related problems appear in [KKS] , 
[KKMS] and [KK]. In the last of these, N. Korevaar and R. Kusner conjecture that 
there is a good moduli space theory for these surfaces. Our methods may be adapted 
immediately to this setting. If M C is an embedded, complete CMC surface with 
A; > 3 ends which satisfies a hypothesis analogous to (5.2) below, then our results im- 
ply that the space of all nearby surfaces of this type, up to rigid motion in M?, is a 
(S/c — 6)-dimensional real analytic manifold. In order to understand the structure of the 
moduli space near surfaces where this hypothesis is not satisfied an argument different 
from the one employed here is needed since, for example, we do not have an analogue 
of the constructions in §6. Recently such an arguement was derived based on the linear 
analysis developed in this paper. Thus we can also assert that the moduli space of such 
constant mean curvature surfaces is locally a real analytic variety, as was claimed in 
[KK]. This arguemnt also applies to give a new, direct proof of Theorem 1.4. These 
results will appear in a forthcoming paper with R. Kusner. 

The outline of this paper is as follows. In §2 we analyze in detail the special case 
when k = 2. The solutions here will be called Delaunay solutions, in analogy with a 
similar family of complete constant mean curvature surfaces in discovered in 1841 
by C. Delaunay [D]. Only in this case may the moduli space be determined completely, 
since, using the symmetry argument discussed above, the equation now reduces to an 
ODE. We also analyze the spectral theory of the linearized scalar curvature operator 
completely in this simple case. §3 collects a number of disparate results about solutions 
of (1.3) with isolated singularities which are used throughout the rest of the paper. This 
contains an explanation of the results of Caffarelli-Gidas-Spruck and Aviles-Korevaar- 
Schoen which state that the Delaunay solutions are good models for arbitrary solutions 
of (1.3) with isolated singularities. We also discuss here the Pohozaev invariants and 
compactness results for solutions g e Ma- The linearization L of the scalar curvature 
operator is studied in §4. We prove Fredholm results for this operator, using and 
expanding earlier work of C. Taubes [T]. Other results here include more detailed 
information on asymptotics of solutions of Lw — 0, as well as the computation, using 
a relative index theorem, of the dimension of the 'bounded nullspace' of L. §5 uses 
these results to establish the structure of AIa near its smooth points. In §6 we show 
that A^A is a real analytic set by writing it as a slice of an 'urmoduli space' Ma with 
the conformal class determined by any g E AiA- Our study of the urmoduli space 
Ma draws on work of A. Fischer and J. Marsden, [FMl] and [FM2]. We also prove a 
generic slice theorem here which shows that slices of Ma by generic nearby conformal 
classes are smooth, even if M.a is not. In §7, our concluding remarks, we discuss three 
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aspects about the nature of M.a concerning which we have not yet obtained satisfactory- 
results. The first of these is the nonexistence of eigenvahies for the hnearization, 
especially for the solutions constructed in [S2] . Secondly, we give a description of local 
coordinates on A^a near smooth points. This is done both on the linear level and 
geometrically in terms of the Pohozaev invariants introduced in §3. At this point we 
also provide a brief discussion of the recent construction [MPU] of dipole solutions for 
the problem. Finally, we address certain natural questions concerning the boundary of 
the geometric compactification M.\. 

The authors wish to thank Rick Schoen for his continued interest and substantial 
advice throughout the course of this work. The first two authors also had a number of 
enlightening conversations with Nick Korevaar, Rob Kusner and Tom Mrowka. 

II. Delaunay Solutions 

In this section we discuss the Delaunay family of solutions. These constitute, up to 
conformal equivalence, the totality of solutions when A has only two elements. 

ODE ANALYSIS 

When k = 2, the condition that A be balanced means that p2 = —pi- It is not 
difficult to show using the Alexandrov refiection argument [CGS] that any positive 
solution of the PDE (1.3) is only a function of the geodesic distance from either pi 
or p2. This equation reduces to an ODE, which takes the simplest form when written 
relative to the background metric dt^ + d6^ on the cylinder M x S""-*^ with coordinates 
{t,0), which is conformally equivalent to S"^\A. Thus, since the cylinder has scalar 
curvature (n — l)(n — 2), and g = u^/^'^~'^\dt^ + dO'^) has scalar curvature n{n — 1), 
u = u(t) satisfies 

d'^ (n-2)2 nin-2) n+2 ^ 

This is easily transformed into a first order Hamiltonian system: setting v — u' (primes 
denoting differentiation by t) 

u — V 

(2-2) (n-2)2 nin-2) n+2 

V = —u ^ ^^"-2 

4 4 

The corresponding Hamiltonian energy function is 

/^ox X (n-2)2 2„ (n-2)2 „ 

(2.3) H{u, v) = -+ ^—^u—- - ^—^u'. 

The orbits of (2.2) remain within level sets of H, and since these level sets are one 
dimensional, this determines these orbits (but not their parametrizations) explicitly. 
The equilibrium points for this flow are at (0, 0) and {u, 0), where 



There is a special homoclinic orbit {uo{t) , vo{t)) corresponding to the level set H = 0; 
it limits on the origin as t tends to either ±00, and encloses a bounded set O in the 
right half-plane which is symmetric across the w-axis, given by {H < 0}. Somewhat 
fortuitously we may calculate explicitly that 

uo{t) = (cosht)^. 

Of course, {H = 0} decomposes into two orbits: this one and the stationary orbit 
(0, 0). It is simple to check that orbits not enclosed by this level set, i.e. those on which 
H > 0, must pass across the v-axis and into the region where u < 0. Thus, since we are 
only interested in solutions of (2.1) which remain positive and exist for all t, it suffices 
to consider only those orbits in Q. Notice that the second equilibrium point {u, 0) is in 
this region, and that all other orbits are closed curves. These correspond to periodic 
orbits {ue{t),Ve{t)), with period T(e), < e < 1. The parameter e may be taken as 
the smaller of the two u values where the orbit intersects the li-axis, so that < e < u 
(note that, strictly speaking, the orbit with e = corresponds to the equilibrium point 
at the origin, but by convention we set it equal to the one previously defined). This 
ODE analysis is also described in [S3] . 

The corresponding metrics on R x S"^"^ (or S'^\{pi,P2}) have a discrete group of 
isometrics, given on the cylinder by the translations t 1— > t + T(e). They interpolate 
between the cylindrical metric u^^^^~'^\dt^ + dO'^) (which is rescaled by the power of u 
so that its scalar curvature is n(n — 1)) and the solution corresponding to the conformal 
factor tto(t). This later solution is nothing other than the standard round metric on 
S"\{pi,p2}j which is therefore incomplete and not, strictly speaking, in the moduli 
space M.A- In all that follows we shall adopt the notation 

(2.5) = uf^ {df + ^6*2) 

when referring to these Delaunay metrics. 

Of greatest concern is the Laplacian for the metrics ^f^; in terms of the coordinates 
(t, 9) on the cylinder we may write this operator as 

4 ^ (n+2) 4 

We use dt for the partial derivative with respect to t, etc. Also, Ag is the Laplacian on 
the sphere S"^^^ of curvature 1. We use the convention that —A is a positive operator. 
Also, we write u = throughout this section. 

It is convenient to replace the variable t by a new variable r(t), depending on e, 
which represents geodesic distance with respect to along the Vt integral curves, 
which are already geodesies (with varying parametrization) for each of the metrics ge- 

2 

r(t) is defined by setting dr/dt — and r = when t — 0. Then 
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(2.6) 



and consequently, using the first equality in (2.6), 

2 — 2ti / 271 — 2 \ 4 

= U '^-'2 dr [U ^-'-^ drj + U ^-^Aq 

^2 2n-2drU^ ^ 

= dl + -—dr + u r.-2Ae. 

n — 2 u 

The function u is still periodic with respect to r, but the period i?(e) now is better 
behaved than the period T(e) above. In fact, whereas 

27]- 

lim T(e) = oo, and lim T(e) = , 

e^O e^u — 2 

(the latter equality is proved by linearizing (2.2) at (•u,0)), we have instead 

27r 

(2.8) lim i?(e) = TT, and lim i?(e) = 



SO that the period stays within a compact interval in the positive real axis as e varies 
between its limits. 

It is somewhat amusing that when n = 4 we can solve (2.1) explicitly, once the 
transformation from t to r is effected; this was pointed out to us by R. Schoen. In fact, 
now di = udr and (2.1) becomes 

(2.9) udr{udru) -u + 2u^ = 0. 

Setting V = udrU as the new dependent variable and u as the independent variable we 
get that vduV — u + 2u^ = 0. Integrating to solve for v, and then integrating again 
leads to a general expression for u. Adjusting the constants we finally get that 



(2.10) „^(^)^^i + (l_e2)cos(2r). 

The parametrization has been arranged so that ^(e) attains its minimum when r = 
7r/2 + Itt, £ G Z. Notice also that this shows that the period R{e) = n when n = 4. 

It is perhaps also instructional, although not necessary within the context of this 
paper, to consider the space of solutions to the problem on the real projective space 
with one point removed, ]RP"\{p}. We can, identify any solution on this space with 
a solution on S"'\{p, —p} invariant under the reflection which exchanges these two 
antipodal points. Transforming the twice-punctured sphere to the cylinder, we are 
looking for solutions u such that u{t, 6) = u{~t, —0). Any such solution is, of course, 
independent of 9 and in the Delaunay family. For each value of the Delaunay parameter 
e there are two possibilities, one taking its minimum value at t = and the other taking 
its maximum value there. These are connected via the cylindrical solution by letting e 
increase to u. The moduli space is thus a copy of M, with the cylindrical solution at the 
'origin.' Alternatively, this moduli space is just the part of the tt-axis in the set H < 
in the {u, f )-plane. As e tends to zero on one end of this moduli space, the solution 
tends to the round spherical metric. As e tends to zero at the other end, the solution 
tends to zero. 

7 



Spectral Theory for the Delaunay Solutions 



Although the Alexandrov reflection argument shows that we have aheady described 
the fuU moduh space M.k when A has only two elements, we proceed further here to 
analyze the linearization of the scalar curvature operator around a Delaunay solution 
g^. This case will serve as the model, and an important ingredient, for the more general 
linear analysis later. Thus, for v a suitably small function, let 

Ne{v) = A,{1 + v)- ^ ^ \ l + v)+ ^ ^ ^ 1 + v)^ 
= A^v + nv + Q{v) 

where Q{v) is the nonlinear, quadratically vanishing term 

/^-.^x ^/ X n(n — 2) ( . . n+2 n + 2 

(2.12) Qiy) = ^ ^ ' -—y 

Solutions of N^{v) — correspond to other complete metrics on the cylinder with scalar 
curvature n(n — 1), and hence correspond to other Delaunay solutions. The linearization 
of about f = is thus given by 



(2.13) ^^^^^^5 



CT=0 



Frequently we shall omit the e-subscript from Lg when the context is clear. 

Using (2.6) and (2.13) and introducing an eigendecomposition {il^j^Xj} for Aq on 
§"-1 we decompose L into a direct sum of ordinary differential operators Lj with 
periodic coefficients on M. The spectral analysis of each Lj, hence of L itself, is ac- 
complished using Bloch wave theory, cf. [RS]. One conclusion is that spec(L) is purely 
absolutely continuous, with no singular continuous or point spectrum. Moreover, each 
Lj has spectrum arranged into bands, typically separated by gaps; L itself has spectrum 
which is the union of all of these band structures: 

spec(L) = [Jspec(Lj). 



Prom (2.13) it is clear that spec(— L) is bounded below by — n. We proceed to analyze 
this spectrum near 0. 

The first question is to understand the Jacobi fields, i.e. elements of the nuUspace of 
L. Any solution of Lij) = may be decomposed into its eigencomponents 

where the aj solve Lj{aj) = 0. Some of these solutions, although not necessarily all, 
may be obtained as derivatives of one-parameter families of solutions of N{v) = 0. It 



is common in geometric problems, cf. [KKS], for the Jacobi fields corresponding to low 
values of j to have explicit geometric interpretations as derivatives of special families 
of solutions of the nonlinear equation. This is the case here, at least for the eigenvalues 
Ao = and Ai = • • • = = n — 1. 

For j = we look for families of solutions of N{v) = which are independent 
of the §^-^ factor. There arc two obvious examples, one corresponding to infinitesi- 
mal translations in t and the other corresponding to infinitesimal change of Delaunay 
parameter: 

e; rj) and 
Now define 

(2.15) 0,=0,-, = ^$,(t,e; 77)1^^0, j = 1,2, 

so that = 0, j — 1,2. Of course, neither (pi nor (j)2 are in since L has no point 
spectrum. In fact, differentiating Ue{t + T{e)) = Ue{t), first with respect to t and then 
with respect to e, shows that 

(2.16) Mt + Tie)) = Mt), Mt + Tie)) + 0i(t)T'(e) = Mt)- 

The first of these equalities states that (pi is periodic while the second shows that 
02 increases linearly, at least so long as T'(e) ^ 0. Hence the (pj are generalized 
eigenfunctions for L with eigenvalue A = 0; their existence and slow growth imply that 
is in the essential spectrum of Lq, hence that of L. It is also not hard to see that 
(pi and (p2 are linearly independent when e ^ u. In fact, if we translate so that it 
attains a local maximum at t = 0, then (pi{0) = 0, whereas (p2{0) = {de/de)/e = 1/e, 
so that these functions are not multiples of one another. When e — u, (pi = 0, but (p2 
does not vanish since 02(0) = 1/u. A second solution of L(p = in this case is obtained 
by translating (p2 by any noninteger multiple of its period. We shall prove below that 
these two functions form a basis, for each e, for all temperate solutions of Lcp = on the 
cylinder; any other solution of this equation must grow exponentially in one direction 
or the other. 

To find Jacobi fields corresponding to the first nonzero eigenvalue of Aq, it is easiest 
first to transform the background cylindrical metric to the round spherical metric. The 
conformal invariance of the equation means that one-parameter families of solutions 
may be obtained from one-parameter families of conformal maps of E>^. For example, 
$1 in (2.14) already corresponds to the composition of g^, with the family of conformal 
dilations fixing the two singular points. We can also consider composition with a family 
of parabolic conformal transformations which fixes one or the other of the singular 
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(2.14) 



Ue{t + ri) 

Ueit) 
_^ Ue+nit) 
Ueit) 



points. Derivatives of these families lead to the new Jacobi fields. To compute them 
we first transform the equation yet again so that the background metric is the flat 
Euclidean metric on M"^ and the singular points are at and oo. The equation becomes 

(2.17) Amr^w + -^-^ — -wr^-^ = 0, 

where the function w on R"^ is related to u on the cylinder by the transformation 

(2.18) v{t,0) = e^''-^^'/^w{e-'), w{p,0) = p^^'""^/^ v{- log p, 9), 

in terms of the polar coordinates (p, 6) on M.^. The Delaunay solutions correspond 
to functions = w^{p). A parabolic transformation of S"^ fixing one of the singular 
points corresponds to a translation of if that singular point corresponds to the 
point at infinity. Infinitesimal translation is just differentiation by one of the Euclidean 
coordinates xj. But 

^0 

dxj ^ dp ' 

and 9j = Xj/p is an eigenfunction for with eigenvalue n — 1. So dw^{p)/dp is a 
solution of Li(p = 0, where Li is written in terms of p instead of t. Transforming back 
to the cylinder once again we arrive at the solution 

, , , , fu'(t) n-2 

(2.19) ) ^ ' 

= e* (0i(t) + ^ 

To obtain the other, exponentially decreasing, solution to Lxil) — 0, we simply observe 
that inversion about the unit sphere in corresponds to reflection in t about t = 0. 
If tig is positioned so that it assumes a maximum at t = 0, then itg(— t) = uj(t)^ and 
we obtain from (2.19) the solution 

(2.20) 0,(t)^e-*('^-0i(t) 



The solutions </>i . . . 04 are also recorded and employed in [AKS] . 

It is also important to have some understanding of the continuous spectrum near 0; 
we do this for each Lj individually. 

Lemma 2.21. For every e G (0, ?i] and j > 1, ^ spec{Lj). Hence the spectral 
analysis of L near the eigenvalue reduces to that of Lq. The only temperate solutions 
of LqcP = are precisely the linear combinations of the solutions (pi and (p2 of (2.15). 

Proof. The eigenvalues of Aq are Aq = 0, Ai = 1 — n, etc. For j > n (numbering the 
eigenvalues with multiplicity, as usual) Aj < —n — 1. Thus 

Lj = u "-2 + 2u "-2 utOt + u "-2 Xj + n 

(2.22) 4 o n+2 4 4 

<u "-2 9^ + 2tt "--^ utOt + n{l — u "-2)-tt "-2 
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Since < tt < 1 for every e and t, the term of order zero in this expression is always 
strictly negative. Standard ODE comparison theory now implies that an arbitrary 
solution of Lj(j) = must grow exponentially either as t — > +oo or — oo. Clearly then 
^ spec(Lj) for j > 1. 

The same conclusion is somewhat more subtle when j = 1. First of all, observe that 
because (f)^ and 04 above span all solutions to the ODE Liip = 0, and since both of 
these functions grow exponentially in one direction or the other, it is evident that is 
not in the spectrum of Li. We wish to know, though, that is actually below all of 
the spectrum of —Li. This is no longer obvious since the term of order zero in 

(2.23) Li = u~^dt + 2u~^utdt + (1 - n)u~^ + n 

is no longer strictly negative. Conjugate Li by the function it^, where p is to be chosen. 
A simple calculation gives 

u-^LivF = u-^dl + {2p + 2)uutdt + A, 



where 



Somewhat miraculously, upon setting p= 2/(n — 2) this expression reduces to 

2n / o {n-2f o (n - 2)^ 2n \ 

(n-2)2 \ * 4 4 J 



where H{e) is the Hamiltonian energy (2.3) of the solution u^{t). This is a negative 
constant for every e G (0, u], and so once again we have obtained an ordinary differential 
operator with strictly negative term of order zero, u 

^2/(n 2) unitarily 

equivalent to Li, and by the previous argument, the spectrum of the first of these 
operators is strictly negative, and in particular does not contain 0. Hence the same is 
true for Li and the proof is complete. 

We can say slightly more about the spectra of the Lj. In fact, when e = tt. 

Til 

and 



spec(— L,) 



— n ( 1 H ) , oo 
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In particular, spec(— Lq) = [— n, oo) and spec(— Li) = [n/ (n — 2), oo). The spectrum of 
—Lu is the union of these infinite rays. As e decreases, gaps appear in these rays, and 
bands start to form (in fact, by (2.8) the first gap appears at — 3n/4). As e decreases, 
the first band extends from — n to some value to the left of — 3n/4; the second band 
begins somewhere to the right of this point, and always ends at 0. It is not important 
for our later work that is always on the end of the second band (or at least, at a point 
where the second Bloch band function has a turning point), but it is rather amusing 
that we may determine this. The explanation is that from (2.16), whenever T'(e) 7^ 
(which certainly happens for almost every e < u), the Jacobi field 4)2 grows linearly. 
This signifies that is at the end of a band, or at the very least, at a turning point 
for one of the Bloch band functions for Lq which parametrize the bands of continuous 
spectrum. As e continues to decrease to these bands shrink further and further. It 
is more revealing to consider the operator written in the geodesic coordinates (r, 9) 
as in (2.7). As noted earlier, as e — > this operator remains periodic, but develops 
singularities every tt units. Geometrically, the metrics g^^ are converging to an infinite 
bead of spheres of fixed scalar curvature n{n — l). The spectrum of —L for this limiting 
metric is a countable union of the spectrum of —A — n on one of these spheres, i.e. 
consists of a countable number of isolated points {— ?7,, 0, n + 2, . . . }, each with inThe 
bands of continuous spectrum have coalesced into these infinite-multiplicity eigenvalues. 
Note that in this limit, the infimum of spec(Li) has increased to zero at a rate which 
may be estimated by a power of e. 

III. Generalities on Solutions with Isolated Singularities 

In this section we collect various results concerning solutions of the singular Yam- 
abe problem with isolated singularities, particularly in the context of conformally fiat 
metrics, which will be required later. 



The Delaunay solutions discussed in the last section are interesting explicit solutions 
in their own right. However, their importance is due to the fact that they are asymp- 
totic models for arbitrary solutions of the singular Yamabe problem, at least in the 
conformally flat setting. Before stating this result more carefully we introduce some 
notation. The half-cylinder (0, oo)^ x S^""*^ is conformally equivalent to the punctured 
ball ]B"^\{0}; an explicit conformal map is given by sending {t,9) to the polar coor- 
dinates {p,0), where p = e~*. A Delaunay solution Ue{t) may be transformed to the 
function 



Using the conformal invariance of the conformal Laplacian, this function solves 



ASYMPTOTICS 



(3.1) 



Ue{p) = P 2 We(-logp). 



(3.2) 
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Theorem 3.3. Let u e C°°(B"'\{0}) he a positive solution of equation (3.2) with a 
nonremovable singularity at the origin. Then there exists some e G (0, u\ and A> Q so 
that u is asymptotic to the modified Delaunay solution in the sense that 

(3.4) zi(p,^) = (l + 0(p°))S,(Ap,^), 

for some a > 0; A corresponds to a translation of the t parameter. The analogous 
estimate still holds whenever u and Ue o,fG replaced by any of their derivatives {pdpYdQ. 

This result was proved by Caffarelli, Gidas and Spruck [CGS] using a fairly general 
and complicated form of the Alexandrov reflection argument. They actually only give 
a somewhat weaker estimate where 0(/9") is replaced by o(l). An alternate argument, 
relying on more direct geometric and barrier methods, was obtained at around the 
same time by Aviles, Korevaar and Schoen in the unpublished work [AKS], and they 
obtained the stronger form. This last work is very close in method to the proof of 
the analogous result for complete constant mean curvature surfaces embedded in R*^ 
presented in [KKS] . It is possible to calculate the decay rate a in terms of spectral data 
of the operator L^, and we shall indicate this argument in the next section. To our 
knowledge it is unknown whether some form of this result has been proved when the 
background metric is not conformally flat, but it seems likely to be true. 

Note that the estimate (3.4) may be restated for the transformed function tt(t, 9) = 
e-("-2)t/2 on the cylinder. Now 

(3.5) w(t,^) = (l + 0(e-"*))w,(t). 

This is the form we shall use. It states that an arbitrary solution on the half-cylinder 
converges exponentially to a Delaunay solution. 

PoHOZAEV Invariants 

A key ingredient in Schoen's construction [S2] of solutions of (1.3) with isolated sin- 
gularities is his use of balancing conditions for approximate solutions. These conditions 
follow from the general Pohozaev identity proved in [S2] . In the present setting, when 
(Q, g) is a compact, conformally flat manifold with boundary, T is the trace-free Ricci 
tensor of ^, X is a conformal Killing field on f], and R{g) is constant, then 

(3.6) / T{X,v)da = {). 

Here u is the outward unit normal to dQ and da is surface measure along this boundary. 
There is a more general formula involving an integral over the interior of Q when the 
scalar curvature is not constant, or when X is not conformal Killing. 

When f2 §"\{pi, . . . there are many conformal Killing fields to use. An 

important class of these are the 'centered dilations.' Any such X is equal to the gradient 
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of the restriction of a linear function £{q) = {q,v), where v e R'^^^ and §" C W^'^^ is 
the standard embedding. Exphcitly, 

(3.7) Xg = v-{q,v)q, 9 G 

since the unit normal to TqS^ is just q. Let c denote the set of all such X. It is a 
subspace of the Lie algebra, o{n + 1, 1), of the conformal group 0{n + 1, 1). 

By (3.6), each X G o(n + 1, 1) determines an element of H^~^{Ei^\A,'R.), associating 
to a hypersurface E the number 

(3.8) ViJ:,X;g)= [ T{X,iy)da 

The dual homology space is generated by the classes of hyperspheres Ej, where, for 
each Pi e A, Sj = dBr{pi) for r sufficiently small, is chosen so that no other pj 
is in the same component of S'^\Sj as Pi. These classes satisfy the single relation 
[El] + ■ ■ ■ + [Sfc] = 0. The number in (3.8) wiU be written Vi{X; g) (or simply Vi{X)) 
when Ti = Tii. Alternately, as suggested by this notation, we shall also regard each Vi 
as a linear functional on the X' i.e. as an element of o*(n + 1, 1). These functionals 
will be called the Pohozaev invariants of the solution metric g. They satisfy 

(3.9) Vi + --- + Vk=^. 

Theorem 3.3 makes it possible to calculate at least some components of the Vi, 
by computing the Pohozaev invariants for the Delaunay solutions; we do this now. 
Since H'^~^ {W^\{pi, — pi}, M) = R, it suffices to compute the single number Vi{X) for 
each X e c. As described above, X may be identified with an element v e R"+^. A 
straightforward calculation, similar to one given in [P2], shows that this invariant is, up 
to a constant, simply the Hamiltonian energy H{e) of the particular Delaunay solution: 

Lemma 3.10. Using the identification of X G c with v G R"'"'"-'^, the Pohozaev invariant 
for the Delaunay metric g^^ on S"\{pi, —pi} is given by 

Vi{X) = CnH{e) {v,p,). 

Here Cn is a nonvanishing dimensional constant (identified explicitly in [P2]). 

The asymptotics in Theorem 3.3 (even with the sharp estimate of a we will give later) 
do not give enough information for the invariant Vi to be computed for every X: the 
point is that if X has associated vector v perpendicular to Pi then Lemma 3.10 shows 
that there is no 'formal' asymptotic contribution to the invariant, but unfortunately 
the decay is not sufficient for there to be no 'perturbation' contribution. The one case 
where this is not an issue is when X is chosen to have associated vector exactly pi. The 
associated invariant Vi{X) in this case will be called the dilational Pohozaev invariant 
and denoted Vi (or Vi{g)) 

Corollary 3.11 ([P2]). The dilational Pohozaev invariant 'Di{g) associated to the 
puncture Pi and the solution metric g is equal to CnH[ei), where is the Delaunay 
parameter giving the asymptotic model for the metric g atpi, as provided by Theorem 
3.3. 
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Local Compactness Properties of the Moduli Space 



It is natural to determine the possible ways that solution metrics g of our problem 
can degenerate. Phrased more geometrically, we wish to determine what the ends of the 
moduli space A^a look like, and to determine a geometrically natural compactification. 
Any A^A, when A contains just two elements, may be identified with any other, and 
we may call this space simply J\42- In this space, all g G are Delaunay, and indeed 
M.2 may be identified with the open set fl C M?, Q = {H < 0} as described in §2. This 
space is two dimensional, smooth, and locally compact. A sequence of elements in M.2 
degenerate only when the Delaunay parameter e (or 'neck-size' of the solutions) tends 
to zero. Recalling that the Hamiltonian energy H is strictly monotone decreasing for 
e G (0, u], this may be restated as saying that Delaunay solutions degenerate only when 
H{e) 0. 

A similar statement is true when A has more than two elements. As we have seen 
in Corollary 3.11, the Pohozaev invariants Vi, and indeed just the dilational Pohozaev 
invariants Vi, determine the Hamiltonian energies of the model Delaunay solutions for 
the metric g at each puncture Pi. Clearly, then, if there exists some particular pi such 
that for a sequence of solution metrics gj G A4\ the dilational Pohozaev invariant 
Vi{gj) tends to zero as j — > 00, then this sequence should be regarded as divergent 
in A^A- The complementary statement, generalizing the situation for AI2) is also still 
valid, and was proved by Pollack [P2]: 

Proposition 3.12. Let gj be a sequence of metrics in Ma, such that for each i = 
1, . . . ,k the dilational Pohozaev invariants Vi[gj) are hounded away from zero. Then 
there is a subsequence of the gj converging to a metric g G Ai\. The convergence is 
uniform in the C°° topology relative to g, or indeed relative to any of the gj, on compact 
subsets o/S^\A. 

The question of what happens to divergent sequences in A^a may also be deter- 
mined. Before describing this we return to Ai2- There is an 'obvious' compactification, 
/A2, which is the closure of O in R^. M.2\M.2 decomposes into two disjoint sets: the 
first contains the single point {0, 0} G M^, which corresponds to the conformal factor 0, 
while the second is the orbit {uo{t) , vo{t)} passing through {1,0}, which corresponds 
to the incomplete spherical metric on S"^\{pi, —pi}- These latter points on the com- 
pactification of Ai2 may be identified with the nonsingular round metric on itself, 
and these are themselves metrics of constant scalar curvature n(n — 1) on S"^. M.2 
is a stratified space; its principal stratum is M2 itself, the codimension one (and one 
dimensional) stratum consists of copies of the nonsingular round metric on and 
finally the codimension two (zero dimensional) stratum consists of the single trivial 
solution of the PDE, which is a completely degenerate metric. 

There is a similar compactification of At a and corresponding decomposition of At a 
when the cardinality of A is larger than two (and still finite). The following is a corollary 
of the proof of Proposition 3.12: 

Corollary 3.13. Let gj G AIa be a sequence such that for some nonempty subset 
of points Pii, ■ ■ ■ ,Pi^, the invariants 'Dij^{gj), . . . ,Vi^{gj) are uniformly bounded above 
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by some —770 < for all j, and all other T>i{gj) tend to zero. Then this sequence 

has a convergent subsequence converging to a metric g G Ma', i-^- 9 is still singular 
at the points A' = {pj^, . . . ,pj^} , but extends smoothly across the points in A\A' and 
R{g) — n{n — 1). If all 'D{gi,pj) tend to zero as i —>■ 00 then either gi tends to zero, 
uniformly on compact subsets of 5" \ A, or else gi converges to the round metric on 



This means that the compactification Ma contains copies of A4a' for certain subsets 
A' C A. In addition it may contain a set whose points correspond to copies of the 
round metric on S"^, and finally a set whose points correspond to the 'zero metric'. 
Later in this paper we shall give a somewhat better description of Ma once we have 
determined the structure of Ma itself better. Notice also that it must be somewhat 
subtle to determine precisely which subsets A' C A have Ma' occurring in A^a- The 
subsets A' which can arise in the description above are determined by the Pohozaev 
balancing condition (3.9). For example, when A has two elements, there is no piece of 
the boundary corresponding to Ma' for A' having just one element; this is because no 
complete solution of our problem exists on S"^\{p}. From the present perspective, this 
holds because T>i{gj) — T>2{gj) = for any j (here A = {^1,^2}) by virtue of (3.9), so 
that if one of these numbers tends to zero, the other must as well. (The difference in 
the signs here from (3.9) is because in (3.9) the same conformal Killing field X is used 
in each Pohozaev invariant, whereas here we use X corresponding to pi for one and 
—X corresponding to p2 = —pi for the other.) The subtle point is that the complete 
Pohozaev invariants, rather than just the dilational ones, are required when A has more 
than two points, and these are not determined just linear algebraically by the location 
of the Pi and the Delaunay models at these punctures. 

IV. Linear Analysis on Manifolds with Asymptotically Periodic Ends 

In this section we prove various results concerning the analysis of the Laplacian 
and the linearization of the scalar curvature operator about solution metrics g E Ma- 
Many of the basic results hold more generally, e.g. for the Laplacian on manifolds with 
asymptotically periodic ends. 

The Linearized Operator 

For the remainder of this paper, L will always denote the linearization about the 
constant solution v — of the nonlinear operator 



(4.1) 



N{v) 




(1+^;) + 



n(n — 2) 



4 



(1+^) 



n + 2 
n-2 



AgV + nv + Q{v), 



where Q{v) is the same as in (2.12), so that, as before. 



(4.2) 



L = Ag + n. 
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Let Ej denote a neighborhood of the puncture pj which is conformally equivalent 
to a half-cyhnder [0,oo)t x We fix these cyhndrical coordinates around each 

Pj. By virtue of Theorem 3.3, L can be treated on each Ej as an exponentially small 
perturbation of the corresponding operator L<=^. for the periodic Delaunay metric g^^ 
which is the asymptotic model for g on Ej. Thus on each Ej we may write 

(4.3) L = . + e-«*F, 

where F is a second order operator with coefficients bounded in C°° as functions of 

it,e). 

The linear analysis of the Laplacian on manifolds with asymptotically periodic ends 
is remarkably similar to that for manifolds with asymptotically cylindrical ends, as 
detailed for example in [Me]. In particular, the Fredholm theory for such an operator 
on exponentially weighted Sobolev (or Holder or ...) spaces has an almost identical 
statement in either case, although the proofs are rather different. The Fredholm theory 
in this asymptotically periodic setting was previously developed by Taubes in [T], 
although we proceed somewhat further into the linear analysis as we need more detailed 
results. It is also possible to develop a full scattering theory for the Laplacian on these 
manifolds. We shall actually require and develop some scattering theoretic results to 
clarify the nature of the moduli space M.a- 

The Fourier-Laplace Transform 

The basic tool for the parametrix construction, upon which all the linear analysis 
relies, is the Fourier-Laplace transform, in a form employed already by Taubes [T] . We 
proceed to develop some properties of this transform, here denoted J^. 

The function spaces we shall use here are exponentially weighted Sobolev spaces 
based on L^iS'^XA; dVg); these will be written iy*(§^\A), or just H?,, for 7, s e M, 
s > n/2. The last condition ensures that the spaces behave well under nonlinear 
operations. To define them, decompose S"^\A into the union of the ends Ei, . . . , Ek 
and a compact piece K. Over K an element h e restricts to an ordinary 

function. Over Ej, h = e^^h, where h e if"([0, 00) x ^'^-^.dtdO). Note that the 
measure here is uniformly equivalent to the one induced by g for any g & AAk- 

The transform JF is, strictly speaking, defined for functions on the whole cylinder 
C = M X S"^~^. It would be somewhat more natural to first develop its properties 
acting on e.g. the Schwartz space 5, but we shall specialize immediately to functions 
with support on just one half of C. So let h{t, 9) E on Ej., and assume /?. = for 
t < 1. Set 

00 

(4.4) h{t,C,0) = J^{h)= J2 e-'''^h{t + k,d). 

k= — (x> 

Assume for the moment that h is smooth. Then, since h decays like e^*, this series 
converges provided Q( = u < —7. We have set C = A* + so that = fi. h{t, C, 0) 
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depends holomorphically on ( in the region i/ < —7. When h is only assumed to be 
in H^, h{t, (, 9) will still depend holomorphically on C, in the same region, but as a 

fimction with values in the space . h is continuous in z/ < —7 as a function with 
values in . These results follow from the Plancherel formalism. 

The tranform T is invertible, and its inverse is given by contour integration. To 
make the following equation clearer, assume that t e R, and i is its reduction mod 1 
(so that < t < 1). Then when £<t<l+lso that we may write t = i+ £, 

(4.5) h{t, 9) = ^ r e''^h{t, C, 9)df,. 

In this formula we integrate along a line Q'C = uq. By Cauchy's theorem this contour 
may be shifted to allow i^q to be any number less than —7. If, as we are assuming, h 
vanishes for t < 0, then h{i, (, 9) not only extends to the lower half C plane, but decays 
like there. Shifting the contour arbitrarily far down shows that the integral (4.5) 
vanishes for any f < 0, as it should. By a similar argument, if /i(t, 6) is defined by the 
integral (4.5) taken along some contour = z/q, where the integrand /i(t, 9) is only 
assumed to be defined along that line, then h G Htj^^^- In particular, if h is holomorphic 
in some lower half plane v < —7, and continuous with values in as a function of t 
up to this upper boundary, then h e H^. 

Next, reindexing the sum defining h gives 

(4.6) h{t + l,C,9)^e^'^h{tX.9). 

This just means that h{t, (, 9) is a section of the flat bundle on §^ x S"^"^ with holo- 
nomy C, around the §^ loop. This bundle is isomorphic to the flat bundle with trivial 
holonomy; the bundle map is given by conjugating by e*^*. Thus the function 

(4.7) h = e-"^*h{t,C,9)e''^^ 
satisfles h{t + 1, C, 9) = h{t, C, 9). 

Fredholm Theory 

The basic Fredholm result for the linearization L may now be stated and proved. 

Proposition 4.8. There exists a discrete set of numbers F C IR such that the bounded 
operator 

(4.9) L : H'+^{M) H^{M) 

is Fredholm for all values of the weight parameter 7 ^ F. In particular, e F, so the 
map (4-9) is not Fredholm when 7 = 0, i.e. on the ordinary unweighted Sobolev spaces, 
but is Fredholm for all values of 7 sufficiently near, but not equal to zero. 

Proof. It suffices to construct a parametrix G for L so that LG — I and GL — I are 
compact operators on and iy^"*"^, respectively. As usual, G may be constructed 
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separately on each piece of the decomposition of S'^\A into a compact piece and the k 
ends El, . . . Ek around each pj . The parametrix construction on the compact piece is 
the standard microlocal one since L is elliptic. We construct a parametrix on each Ej 
using the Fourier-Laplace transform. 

Fix one of these ends, Ej, and let be the model asymptotic Delaunay metric 
for the fixed metric g there. The corresponding model operator L^, as in (4.3), has 
periodic coefficients (for notational convenience wc assume that the period is one here) 
so it acts on sections of the fiat bundle with holonomy C described above by the obvious 
rule (L^h) = L^h. This induced operator looks just like in local coordinates (t, 6). 
This step is the same as conjugating by JF. We may proceed further and conjugate 
Le by e*^* so as to act on the trivial flat bundle. This final induced operator, which 
depends holomorphically on C, will be called -£'e(C): 

(4.10) Le(C) = e-'^'L,e''^\ 

The main point of the proof is that Lei.0 inverse (say on L^(§-^ x S"'"^)) 

which depends meromorphically on Q. This will be a direct consequence of the analytic 
Fredholm theorem, which is proved for example in [RS] and states that a strongly 
holomorphic family of Fredholm operators, depending on the complex variable either 
fails to be invertible for every C, or else is invertible for all C, except for those in some 
discrete set in the parameter space. To check that this result is applicable, simply note 
that I'e(C) is elliptic for every C, (as already observed) depends holomorphically on C, 
and, by the various transformations we performed, acts on a fixed function space (on a 
fixed bundle) over a compact manifold. Hence it forms a family of Fredholm operators; 
standard considerations show that the holomorphic dependence of the coefficients of 
this operator on C, ensure the strong holomorphy of the family. 

Our next task is to show that L^{C) is invertible for some value of C- Once this 
is accomplished, the analytic Fredholm theorem will imply that the set of poles {Cj} 
of Le(C)~^ is discrete in C. The invertibility of L^{C,) at some C, is equivalent to the 
invcrtibility of acting on the flat bundle with holonomy C,. When C, is real, this 
operator is self-adjoint, hence invertible so long as it has no nuUspace. A solution of 
LeV' = on this bundle lifts to a function ip{t, 6') on R x satisfying ■il){t +1,6) = 
e'^^il){t,9). Since C is real here, this lift is bounded and 'quasi-periodic' However, we 
have already shown in §2 that such a function can exist only when ^ = 0; in fact, 
the arguments there show that any temperate solution must be constant on the cross- 
section S"^~^, i.e. be a solution of the reduced operator Lq. In addition, the functions 
01 and 02 of (2.15) are independent solutions of this operator, hence span the space 
of all solutions since Lq is a second order ODE. Finally, neither of these transform 
by a factor e**^ over a period except for C = 0. Thus is invertible for every real 
C G (0,27r). This proves that L^[C,)~^ exists and depends meromorphically on Q. 

We make some remarks about the set of poles {Cj} of this meromorphic family of 
inverses. First, if C, is in this set, then so is C + 27r£ for any £ e Z. In fact, although 
Le(C+27r) is not equal to -£'e(C)) these two operators are unitarily equivalent: the unitary 
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operator intertwining them is multiplication by the function e . So, for some fixed C, 
either both of these operators are invertible or neither of them are. By discreteness and 
this translation invariance, the set {Q} has at most finitely many inequivalent (mod 
2tt) poles in any horizontal strip a < < b. Note in particular that the only poles 
on the real line occur at the points in 27rZ. By discreteness again, there are no other 
poles in some strip — e < < 0- We shall usually restrict attention then only 

to values of ( in the strip < C < 27r (or equivalently, in C/27rZ); V will denote the 
set of poles in this strip. Finally, note that the adjoint of L^{Q is precisely L^{Q, and 
so V is invariant under conjugation. This means that we can list the elements of V as 
follows: 

(4.11) V = Ve^{..., -C2, -Ci, Co = 0, Ci, C2, . . . }. 
Also define 

(4.12) r = r, = {7, = ^0 : G 7^} = {■ ■ • , -71, 0, 71, • • • }• 

The 7j form a strictly increasing sequence tending to infinity, and in particular, F is a 
discrete set in M. We shall provide an interpretation of the nonzero elements of V in 
the next subsection. 

Next we unwind this inverse L^{(^)~^ = GeiC) obtain a parametrix for on Ej. 
The inverse for acting on the flat bundle with holonomy ^ is given by G = e^^^Ge"*^'' , 
for every C, ^V. To obtain an inverse for we simply conjugate G with the Fourier- 
Laplace transform JF. Recall, though, that this makes sense only if the contour integral 
in the definition of !F~^ avoids V. So, by conjugating, and in doing so, taking this 
integral along a contour = —7, for any 7 ^ F, we obtain an inverse for Lg. Since 
integrating along this contour produces a function in H^{W x S""-^), we have obtained 
an inverse G = G^^^ for acting on functions on Ej, supported in t > 1 for all 
7 ^ F. Clearly this G is a parametrix with compact remainder for L = + e~"*F 
acting on H^{Ej). 

As a final step, couple these parametrices we have produced on each end Ej to the 
interior parametrix to obtain a global parametrix for L on all of S'^yA with compact 
remainder. Both right and left parametrices of this type may be obtained this way, 
and their existence implies immediately that L is Fredholm on whenever 7 ^ F. It 
is a simple exercise to check that L does not even have closed range when 7 e F, in 
particular, when 7 = 0. 

It is obviously of interest to determine when L is actually injective or surjective. In 
general this is a rather complicated question, but we make note of the following result. 

Corollary 4.13. Suppose that L has no Lp' nullspace (i.e. is injective on for^ < Oj, 
then for all S > sufficiently small ( with S < a) 

L : Hp'^ — > Hg is surjective, 

L : — > H^g is injective. 

This first statement follows from the second and duality, since L is self-adjoint on 
L^, i.e. when 7 = 0. 
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Asymptotic Expansions 



The Fourier-Laplace transform may also be used to obtain existence of asymptotic 
expansions for solutions to Lw = on each end Ej. For our purposes in this paper, it 
will only be important to know that any such w has a leading term in its expansion, 
which decays (or grows) at some specified rate, and an error which decays at a faster 
rate. One use of this will be to estimate the exponent a appearing in (3.5). However, 
the full expansion is not much harder to prove, so we will do this too. 

The starting point is the fact, discussed above, that the Fourier-Laplace transform 
of a function w e H^(R x S"'"-'^), supported in t > 0, say, is a function w{t, 9, () which 
is holomorphic in the half-plane Q'C < —7 (and taking values in if*(S^ x S"^"^)). For a 
general function of this type, the most that can be concluded is that it has a limit on 
the line Q'C = ~7 the appropriate L^-sense discussed earlier. 

When the solution w is transformed to a function w{r, 9) — w{— logr, 9) on B"\{0}, 
elliptic regularity shows that is a conormal distribution with respect to the origin, 
i.e. has stable regularity when differentiated arbitrarily often with respect to the vector 
fields rdr, dg (this is just the same as w itself having stable regularity with respect to 
dt, dg); [Ml], for example, contains a discussion of conormal regularity. For degenerate 
operators of a type closely related to L, one expects solutions to be poly homogeneous 
conormal, cf. [Ml] or [Me]. Polyhomogeneity is simply the property of having an 
asymptotic expansion in increasing (possibly complex) powers of r and integral powers 
of logr, with coefficients smooth in 9 (as functions of t these expansions are in analogous 
powers of and t). Alternately, polyhomogeneous conormal distributions may be 
characterized as those with Mellin transforms, already defined and holomorphic in 
some lower half-plane in C, extending meromorphically to the whole complex plane, 
with only finitely many poles in any lower half-plane, all of which are of finite rank. 

w will not be polyhomogeneous, except when e — u and the underlying metric 
is cylindrical. In fact, iv (cut off so as to be supported in r < 1) will have Mellin 
transform still defined and holomorphic in a half-plane and extending meromorphically 
to all of C, but now its poles are arranged along lattices on a countable discrete set 
of horizontal lines. Each pole is still finite rank. This meromorphic structure of the 
Mellin transform of tu is equivalent to the fact that the terms in the expansion for w 
have the form e~^^Wj^ where Wj is periodic in t. The residues at the poles along a fixed 
horizontal line at height /3 will be the Fourier series coefficients of Wj. 

This more general discussion has involved the Mellin transform as r — > (or equiva- 
lently, the standard Fourier transform as t ^ 00). We shall revert now to the Fourier- 
Laplace transform; it is much less flexible than either of the other transforms, since it 
presupposes periodicity, but suffices here for our immediate purposes. 

Suppose that «; G i/^(M x S'^-i) for some 7 ^ F is supported in t > and solves 
= / for some compactly supported smooth function / (hence w G for every 
s by elliptic regularity). For example, we could take which solves L^cj) = 0, and 
let w = x4> where x is a cutoff function having support in t > and equalling one 
for t > 1. Taking transforms we get Lw = /; this function on the right, f{t,9,Qi is 
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obviously smooth in (t, 9) and entire in Applying the inverse G from the last section 
gives w{t, 9, C) = Gf{t, 9, (). The right side of this equation is meromorphic in ( with 
poles at some subset of points in V, hence the same is true for the left side as well. 
Notice that the poles of G in < —7 must be cancelled by zeroes of / since w is a 
priori known to be regular in this half-plane. 

The function w is recovered by inverting integrating along the line = ~7- 
However, by Cauchy's theorem this integration may be taken along any higher contour 
S^C = —7', so long as the interval [7', 7] does not contain any points of F. Thus, for 
any such 7', w G Hy. If the contour is shifted even further, so as to cross a point 

in V, i.e. a pole of G, then the resulting integral along the line $5^ = —7" produces 
a new function v. Since {L^v)'" — L^Gf — f and / is entire, it follows that v also 
solves Lv = f, and so w and v must differ by an element of the nuUspace of Lg. By 
construction v e i^y/, and so we have decomposed w as a sum v + where L^'^ = 0. 
ifj is given by the residues of w at the points Q G V with —7 < '^{—(j) = —jj < —7". 
(We use — instead of Q for notational convenience only.) If there is just one such 
point, say —Q, then it is clear that in fact '0 G H^.^^ for every e > 0. Actually, it is 
not hard to show, using results from §2, that this ij) must grow or decay exactly like a 
polynomial in t times e^^^. 

This process may be continued by moving the contour past more and more poles. 
At each step, we have decomposed w into a sum of solutions of L^'^ = and a term 
which decays at a rate given by the height of the contour. 

Proposition 4.14. If w solves L^w = f for some compactly supported function f on 
(0, cxo) X S"-"-^, with w G for some 7 ^ —T, then, as t ^ 00, wit, 9) ~ ^ 6*), 
where each i/jj solves L^ifjj = and il^j decays like a polynomial in t multiplying e~^^''', 
where 7^ G — F as in (4-12). 

Note that this result shows that the poles Cj, or at least their imaginary parts, 7^-, 
correspond to the precise growth rates of solutions of L^(p = on the cylinder. 

A similar, though slightly more complicated, expansion holds for the elements of the 
nuUspace of the linearization L in (4.2). 

Proposition 4.15. Let w solve Lw ~ 0, at least on some neighborhood of the puncture 
Pj, and lie in for some 7 ^ — F, where F is the set of imaginary parts of poles 

corresponding to the inverse G for the model for L near this puncture. Then w has 
an asymptotic expansion 

00 

j,k=0 

as t ^ (yo. In this sum,, the 'leading terms' i/jj^q are solutions of L^ijjj Q = corre- 
sponding to the poles Qj above = —7. These decay, as before, like a polynomial 
in t multiplying e~^^*. The higher terms '4>j,k, A; > 0, decay like a polynomial in t 
multiplying e~^'^^^^°-'>* . 

22 



Proof. As before, we may assume that w is supported in t > and Lw = / is compactly 
supported and smooth. By (4.3), 

(4.16) L,w ^ -e-'^^Fw + f. 

Now conjugate by the Fourier-Laplace transform, and apply the inverse G of Lg to get 

(4.17) w = G(-e-"*F«;)^+ Gf. 

The term on the left is holomorphic, a priori, in the half-plane < ~7- The second 
term on the right is entire, while the first term on the right is holomorphic in < —7 
and extends meromorphically to the slightly larger half-plane S^C < —7 -|- a. As before, 
w can be recovered by integrating along = —7 — e; if this contour is moved up 
to = —7 + a — e, and if the strip —7 < < —7 -|- a contains no poles, then 
we find that w e H^_^. If there are poles in this strip, then, exactly as in the last 
Proposition, w decomposes into a sum v + il^ with v G H^_^ and L^if; = 0. •0 in turn 
decomposes into a sum of terms i/jj^q corresponding to the various poles in this strip. 
This improved information may be then fed back into (4.16) and (4.17). Now the right 
side of (4.17) extends meromorphically to the half-plane < —■^ + 2a, and the contour 
may be shifted further to get more contributions to the expansion for w. Continuing 
this bootstrapping yields the full expansion. Details are left to the reader. 

The Deficiency Subspace 

Particularly important for us in the application of the implicit function theorem will 
be the pole E V for G (for any value of the Delaunay parameter e). Specifically, 
we will be concerned with those solutions of L^w — and Lw — which are in 
for every 5 > 0, but not in any H^g. Here we are still only concerned with the local 
behavior of these solutions on each end Ej. Their global nature will be discussed later. 

The solutions to L^(j) = in Hg Ht^ are given in (2.15) as the functions (pi = (pi^^ 
and 4'2 = (p2,e- They depend on t but not on 9. The fact that 4'2 grows linearly in t 
indicates that G has a pole of order 2 at = 0. These are the only temperate solutions 
of Le0 = on the whole cylinder; any other solution grows at an exponential rate (with 
possible exponents given by the values of = T). 

If the solution metric g has model Delaunay parameter e on the end Ej, then these 
functions (for that value of e) may be cut off and transplanted onto this semi-cylinder. 
They do not decay, but it is easy to create a sequence of cut-offs Xi4^ (0 — 0i or (^2) 
with the following properties: each Xi is compactly supported, and has support tending 
to infinity as i tends to infinity. The norm of Xi^P equals one for all i, but the 
norm of L{xi(p) tends to zero (L is the linearization, not one of the models Lg here). 
The existence of such sequence is a standard criterion for showing that L does not 
have closed range on L^. In any event, we think of (pi, (p2 as constituting the bounded 
approximate nuUspace for L (here 'bounded' is loosely interpreted to encompass the 
linearly growing 02)- 
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Define now a finear space W generated by the functions </)i^£ . , 02, e , cut off in a fixed 
way so as to be supported in t > and transplanted on each end Ej; the ej are the 
Delaunay parameters for the particular model metrics g^. on Ej. Since there are k 
ends, W is 2/c-dimensional. For reasons that will become clear in the next section, we 
call W the deficiency subspace for L. Clearly W C Hg for every s and every 5 > 0. 

The Linear Decomposition Lemma 

There are two important corollaries of Proposition 4.15 and its proof which we single 
out in this section. The first concerns the behavior of solutions of the inhomogeneous 
equation Lw — f on each end Ej, while the second studies the exact value of the 
exponent a appearing in Theorem 3.3 and then later in (4.3), etc. 

As already pointed out in Corollary 4.13, if L has no global nuUspace, then we 
can find a solution w G H^'^'^ to the equation Lw = f for every / G Hg, whenever 
5 > 0. In particular, this holds whenever / G H'^g. Clearly, whenever / decays at some 
exponential rate like this, we expect the solution w to be somewhat better behaved 
than a general Hg'^'^ function; of course, it is immediate that it is in this space for any 
5 > 0, but we can do even better. This is the subject of what we will call the 

Linear Decomposition Lemma 4.18. Suppose f G H^g for some 5 > sufficiently 
small, and w G Hg~^^ solves Lw = f. Then w G © W, i.e. w may be decomposed 

into a sum v + (f) with v decaying at the same rate as f and cj) in the deficiency subspace 
W. 

Proof. Clearly this question may be localized to each end Ej, and in this localized 
decomposition will be a combination of 0i and 02- The decomposition is achieved by 
exactly the same sort of shift of the contour in the integral defining JF~^ across a pole 
of G. Here the contour is being shifted from = —5 to Q( = +5; the pole crossed is 
the one at C = 0. 

The second corollary deals with the rate at which a general solution of (1.3) on a 
punctured ball, singular at {0}, converges to the radial Delaunay metric. We assume 
the simpler statement of Theorem 3.3 that u decays to at some exponential rate 
a, and use the linear theory to find the optimal rate. First transform the punctured 
ball to the half-cylinder (0, oo) x §"^~^, and assume all functions are defined here. 
Write the solution w as a perturbation (1 + v)tte, so that N^{v) = 0, where is the 
nonlinear operator (2.11). This equation is the same as L^v = —Q{v), where Q is the 
quadratically vanishing function in (2.12). We already know that v G Hta (for all 
s), hence Q{v) G H^2a- -^Y contour-shifting arguments above, this implies that v 
itself decays at this faster rate, at least provided 2q; < 71 , where 71 is the first positive 
element in T. Continue this process until this first pole at Ci has been crossed; the 
conclusion is that v decays exactly like e~'^^^ (possibly multiplied by a polynomial in 
t). We could also bootstrap further and obtain a complete asymptotic expansion for 
V, of the same general form as in Proposition 4.15 above, although more exponents 
occur because of the nonlinearity. (This is a simple form of the argument in [M2] ) . We 
summarize this discussion as 
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Proposition 4.19. The exponent a occuring in Theorem 3.3 governing the rate of 
decay of a general solution u to its model Delaunay solution is equal to the first 
nonzero element 71 in the set T corresponding to G^. The function u admits a complete 
asymptotic expansion into terms of increasingly rapid exponential decay. 

The Bounded Nullspace 

In this last subsection of §4, we finally come to the global behavior of that portion 
of the nullspace of L corresponding to the pole at 2; = 0. This space, B, which we shall 
call the 'bounded nullspace' of L, is defined by 

(4.20) B = {veHI:Lv = ^,vi H'_g}. 

The full nullspace of L on Hg will be the direct sum of B and the nullspace (which, 
by Proposition 4.15 is the same as the nullspace of L in H^g). In particular, when this 
latter space is trivial, B is the full nullspace of L in H^. By the linear decomposition 
lemma 4.18, B is already contained in H^g © W. The purpose of this subsection is to 
determine the dimension of B. 

As usual, an index theorem is the principal tool for calculating dim(i3). Fortunately 
we require only a relative index theorem, which computes the difi^erence between two 
indices in terms of asymptotic data, rather than global data. To be more explicit, for 
any 7 ^ F set 

(4.21) ind (7) = dim ker L\fjs — dim coker L\jjs ■ 

This is obviously independent of s. Since the adjoint of L on Ht^ is L on i^Z^, duality 
implies that 

(4.22) ind (—7) = —ind (7) for every 7 ^ F. 

If 7i and 72 are any two allowable values (i.e. neither are in F), the relative index 
with respect to these two numbers is simply the difi'erence 

rel-ind (71, 72) = ind (71) - ind (72). 

In particular, using duality again, 

rel-ind (5, -b) =2 ind (5) 

(4.23) =2 ^dimker L|^| — dim ker L|^s^ j 

=2dim(B). 

As noted above, this relative index can be shown, on fairly general principles, to 
be computable in terms of asymptotic data for the operator L. Finding a specific and 
computable formula is another matter and, to our knowledge, there is no general result 
of this sort available for asymptotically periodic operators. However, such a result is 
available for for operators associated to asymptotically cylindrical metrics, and we will 
use this instead. Our result is 
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Theorem 4.24. dim(B) = k. 

Proof. By (4.23), it suffices to show that rel-ind (5, —5) = 2k. By the usual stabihty 
properties of the index (hence any relative index) under Fredholm deformations, we 
compute this number by choosing a one-parameter family of Fredholm operators L'^, 
< o" < 1, with = L, and an operator for which this relative index is computable. 
Since L is just Ag + n, we shall choose a one-parameter family of metrics g'^ and define 
L'^ to be Agcr +n. The metrics will agree with g except on each of the ends Ej, where 
we make the following homotopy. First deform g on each end to its model Delaunay 
metric g^, and then deform each through Delaunay metrics to the cylindrical metric 
g^. The metric g^ will agree with the original g on any large fixed compact set, and 
will equal g^ on each end Ej. This metric is now an exact b-metric, in the language of 
[Me] (more prosaically, it has asymptotically - and in this case, exactly - cylindrical 
ends), and the corresponding operator is an elliptic 6-operator. 

Of course, we still need to prove that L"^ is Fredholm on Hg and H^g for every 
< cr < 1, provided 5 > is sufficiently small. In the part of the deformation 
where g is homotoped to its model Delaunay metric g^ on each end, this is obviously 
true. For the remaining part of the deformation, we need to know that the elements 
of {7j e T,j > 1}, which are the weights for which L will not be Fredholm, remain 
bounded away from as e varies between its initial value and u. This is precisely the 
content of Lemma 2.18. 

Now, to apply Melrose's Relative Index Theorem [Me], note that we may as well 
consider the operator — ^^^^L, which on each end takes the form + dg + {n — 2). 
The set V for this operator (called specfe(L^) in [Me]) is {±iy/n — 2, ±7j,i > 1}, where 
7j are all strictly positive, tending to infinity, and obtained in a straightforward manner 
from spec(Z\0). The prescription to calculate the relative index is to first consider the 
'indicial family,' which in this case equals Aq + (n — 2) — (it is obtained by passing to 
the Fourier tranform with respect to t, which carries dt to —iQ. This is a holomorphic 
family of elliptic operators on S""-*^, and its inverse G(C) has poles exactly at the points 
of V. Since we are computing the jump in the index as the weight changes from —S to 
+5, we need to compute the 'degree' of each pole (as defined in [Me]) for each element 
of this set with imaginary part equal to zero. Each of the poles ±V?i — 2 is simple and 
of rank one, and so the degree of each is also equal to one. Finally we need to sum over 
all poles and over all ends Ej, because the preceding discussion is local on each end. 
For each end there are two poles of G with imaginary part zero, each contributing a 
degree of one to the computation, and there are k ends; The sum of all of these is 2k, 
and this is the relative index rel-ind {6, —6). 

V. The Moduli Space: Smooth Points 

In this section we commence the study of the moduli space Ai\ itself. Here we use 
the linear theory developed in the previous section in a straightforward way to study 
neighborhoods of the 'good points' g e Ma, where the associated linearization has no 
nuUspace. In the next two sections we develop ideas to study neighborhoods of the 
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nonsmooth points. 

The implicit function theorem is directly applicable only when the linearized operator 
Lg = Ag + n is surjective on the appropriate function spaces. These function spaces 
should be tangent to a suitable space of metrics conformal to g and with prescribed 
growth conditions on each end of §"'\A. Let g e M.a, and consider all nearby metrics 

4 

conformal to g of the form (1 + v)^^g. The most natural class of such metrics is 

(5.1) Metl"^^ = {(1 +1;)^^, v e Htg, -y > -1}. 

Clearly TgMet'+/ = H'_Y. The operator Ng{v), as defined in (4.1), can be considered 
as a map from this space of metrics to H^g, and as such is obviously C°°. However, it 
is never surjective, since —5 < 0. The obvious alternative, to consider metrics growing 
at the rate +5 is unsuitable because of the nonlinearity. 

The case where this difficulty is easiest to remedy is when Lg has no nuUspace: 

(5.2) ker (Lg) n l2(§"\A, dVg) = {0}. 

This will be our standing hypothesis in the rest of this section. 

We can rephrase this condition by recalling, by virtue of Proposition 4.15 on asymp- 
totics of solutions of Lcj) = 0, that if were in the nuUspace of L, then it would 
decay at the exponential rate e~^'^* on each end, where 71 = 71 (cj) is the first nonzero 
element of V^j on the end Ej, corresponding to the model Delaunay metric g^. there. In 
particular, under the hypothesis (5.2), for 6 G (0,71), the operator L has no nuUspace 
on H^g for any s G M. Now apply Corollary 4.13 and the linear decomposition lemma 
4.18 to conclude that (5.2) is equivalent to the statement that 

L : ifl+2(§-\A, dVg) ®W ^ Jf^d(§"\A, dVg) 

(5 3) 

^ ' ^ is surjective for any < S < min{7i(ej)} and s e M. 

j 

W is the 2/c-dimensional deficiency subspace introduced in the last section. Because 
we work with a nonlinear equation, we assume that s > n/2. 

Theorem 5.4. Suppose g G Ma is a metric such that the hypothesis (5.2) is satisfied. 
Then there exists a 2k -dimensional open manifold W, the elements of which are func- 
tions V on S"\A to he described below, and with a distinguished element vq G W, such 
that the nonlinear map 

Ng : iy^+2(S'^\A; dVg) X W ^ H^si^^X^; dVg) 

defined by Ng{v,il^) = Ng{v -\- where Ng is the operator given by (4--1), restricted to 

a neighborhood of in the first factor, is real analytic. Moreover, the tangent space 
T„qW is identified with the deficiency subspace W , and using this identification, the 
linearization Lg of the map Ng is surjective, as in (5.3). 
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Corollary 5.5. Suppose g e A^a satisfies the hypothesis (5.2). Then there is an open 
set U C AIa containing g, such that U is a k- dimensional real analytic manifold. 

The proof of the Corollary follows directly from the surjectivity statement (5.3), 
by a straightforward application of the real analytic implicit function theorem. The 
dimension k is of course the dimension of the nuUspace of L as a map (5.3). This 
nuUspace is the bounded nuUspace B defined in the last section, and by Theorem 4.24 
its dimension is k. So to conclude the proof of the theorem and the corollary it suffices 
to construct W, and to show that Ng is a real analytic mapping. 

The only difficulty in the construction of W is that the elements of W are not, 
generally, bounded either above or below, so we could not use (1 + t;)4/("^~2) 
conformal factor for most v G W . However, recall from (2.15) that each element 
j = 1, 2, of the bounded nuUspace for the model problem on the cylinder is the tangent 
to a curve rj —>■ ^j{t,e,r]) of actual conformal factors. The actual curve e W we 
would want to use for each such element of W on an end has the form 



The manifold W is constructed by gluing together these local definitions from each end 
Ej in an essentially arbitrary, but smooth, manner. By construction, TqW = W, as 
desired. 

To check the real analyticity of Ng, the main step is to write Ng{v + il^) in such a 
way that it clearly lands in H^g. Hence, on each end write Ng = + e~^*Q, and then 
use a common formula for the remainder in Taylor's theorem to obtain 



Since N^(ip) = 0, every term on the right here is in Hg so long as 5 < e. Furthermore, 
every term is real analytic in {v,ip). Real analyticity of Ng on the interior, away from 
the ends, is even easier. 

The rest of the proof is now standard. 

VI. The Real Analytic Structure of Ma 

In this section we prove that A^a is always an analytic set by representing it as the 
slice of an infinite dimensional real analytic manifold Ma with the conformal class go- 
We go on to prove a generic slice theorem, that even if Ai\ itself is not smooth, generic 
nearby slices of Ma by other conformal classes will be. 



As described above, we wish to regard A^a as the slice by the standard conformal 
class [^o] of the infinite dimensional set Ma consisting of all metrics, not necessarily 





1 



The Urmoduli Space Ma 
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in the standard conformal class on S"^, with constant scalar curvature n{n — 1). For 
obvious technical reasons, we consider only those metrics which satisfy appropriate 
growth and asymptotics conditions near the punctures pj (this will be be elaborated 
on below.) Consideration of this 'big,' or 'ur-,' moduli space (when the underlying 
manifold is compact) was first undertaken by Fischer and Marsden in the early '70's 
([FMl] and [FM2]), motivated in part by concerns in general relativity. They proved 
that it is a smooth Banach manifold (in any one of of a number of standard Banach 
completions), provided a certain overdetermined linear equation has only the trivial 
solution. 

Their precise set-up was to consider Mp, the set of all metrics (of some fixed finite 
regularity) with scalar curvature function p, a fixed function on the underlying manifold 
X . To show that this set is a Banach manifold, it suffices to show that the linearization 
of the scalar curvature map (which assigns to any metric g its scalar curvature function 
R{g)) is surjective at any g G Mp. Since this map R carries metrics to functions, its 
linearization, which we will call ctg here, carries the tangent space at g of the space of 
all metrics, i.e. the symmetric 2-tensors, to scalar functions. A. Lichnerowicz [Li] had 
earlier computed that, for h a symmetric 2-tensor, 

(6.1) ag{h) =66h- Ag{trh)- {rg,h). 

Here S is the divergence operator on tensors, and Vg is the Ricci tensor for g. It is easy 
to check that ag reduces to (a multiple of) the linearization operator Lg we have been 
studying when /i is a multiple of the metric g (and hence is tangent to the conformal 
class [^f]), when p = n{n — 1). It is straightforward to check that the symbol of ag 
is surjective, so it follows that the symbol of its adjoint a* is injective, and then that 
the symbol of ag o a* is an isomorphism. Therefore, this last, fourth order operator 
is elliptic; at least when X is compact, it is Fredholm, and in particular has closed 
range of finite codimension. Since range {ag) D range {ag o ct*), the range of ag itself is 
closed and of finite codimension. (It does, however, always have an infinite dimensional 
kernel.) By virtue of all this, when X is compact, ag is surjective if and only if a* 
is injective, and this («*(/) = 0) is precisely the overdetermined equation referred to 
earlier. 

From (6.1) it follows readily that 

(6.2) «*(/) = Hess,(/) -^Z\,(/) - /r,. 

If a*{f) — 0, then taking the trace we also get that Ag{f) — —{pj {n—l))f; substituting 
this back into (6.2) we obtain finally the condition 

(6.3) a;{f) = =^ Hess,(/) = - f r, - -^g) f. 



Fischer and Marsden were able to show that (6.3) (or equivalently, (6.2)) has only the 
trivial solution / = except possibly in the special cases where p is a nonnegative 
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constant and pj {n — 1) is in the spectrum of —Ag. Thus, except for these cases, Mp is 
a smooth Banach manifold. 

In the case of interest here, the underlying manifold S"^\A is not compact, and 
p/(n — l) = n is always in spec {Ag). We address the former of these concerns first, 
since the Fredholm theory of §4 may be adapted easily to this problem. As before 
iy*(§"\A) will denote the Sobolev space of scalar functions decaying like e'"'* on each 
end, with respect to fixed cylindrical coordinates {t, 6) there and with respect to a 
fixed metric g E AiA (any element of the moduli space can be used as a fixed back- 
ground metric to define these Sobolev spaces). Define also i/:^(S"\A, S^T*(§"^\A)) to 
be the space of symmetric 2-tensors h = e'^*k, with k in the unweighted Sobolev space 
iy«(§'^\A, S2T*(§^\A); dVg). 

Proposition 6.4. For any g G Ma, the operator (6.1) is hounded as a linear map 

ag : iy^+2(§"\A,52/m2(T*(S^\A))) i?^(§"\A). 

It has an infinite dimensional nullspace, and the closure of its range has finite codi- 
mension for a/Z 7 e R. There exists a discrete set V such that for all 7 ^ F', this map 
also has closed range. 

Proof. The boundedness assertion is immediate. The proof of the closed range part is 
very close to the analogous one given in §4 for the scalar operator Lg, so we just provide 
a sketch. As in [FM2] and above, the surjectivity of the symbol of ag guarantees that 
we merely have to show that Ag = ag o a*, which is elliptic, has closed range for all 
7 ^ r'. The operator Ag is asymptotically periodic on each end Ej of S"^\A, and the 
Fourier-Laplace transform may again be employed to construct local parametrices on 
each of these ends. The fact needed to make this work is that the holomorphic family 
of elliptic operators A^{(j)), on the compact manifold x S"^~^, and associated to the 
Delaunay model operator Ag^_, is invcrtiblc at some (p E C. In this case the analytic 
Fredholm theory, and the rest of the construction in §4, proceeds exactly as before. 
Again we prove this for e R, since for such 0, ^e(0) is self-adjoint. Invertibility is 
equivalent to the injectivity of the induced operator a* {(f)), and this in turn is implied 
(by taking traces) by the injectivity of the model operator L^{(p). The invertibility of 
this last operator for all (/) G (0, 2tv) (we are again assuming the period of a^, et al., to 
be one for simplicity) was already proved in §4, and so we are done. 

To proceed with the program of [FM2], we also have to show that a* is injective on 
S'^\A for any g G A4a- If this were true, then the urmoduli space Ma would be a real 
analytic manifold in a neighborhood of every g G TWa, i-e. there would exist a (local) 
thickening of Ma into the real analytic moduli space Ma such that Ma — Ma fl [go], 
as desired. Unfortimatcly this is not quite so simple, since, as already noted, p/{n — 1) 
is always in the spectrum of Ag] we are thus precisely in the case not treated by [FM2]. 

The equation (6.3) is rather similar to the well-known equation of Obata, which also 
appears in the study of the scalar curvature, and actually reduces to Obata's equation 
when X is an Einstein manifold. In analogy with the situation there - in particular, 
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Obata's characterization [0] of {Ei^,go) as the only compact manifold admitting non- 
trivial solutions of his equation - Fischer and Marsden conjectured that (6.3) admits a 
nontrivial solution, at least when X is compact, if and only if X is the standard sphere. 
This is now known to be false: indeed, let X be the product of the circle of length L, 
§^(L), with any Einstein manifold E'^~^ of positive scalar curvature, e.g. S"'"^. Then, 
if the scalar curvature of E is normalized to be {n — — 2), (6.3) implies that / 
depends only on the variable t along the circle, and satisfies the ODE /" + (n — 2)/ = 0. 
Thus, (6.3) admits a nontrivial solution whenever L is an integer multiple of — 2. 
There are also noncompact counterexamples to their conjecture: the simplest are just 
the universal covers M x E'^~^ of the compact examples above. 

The characterization even of all compact manifolds admitting solutions to (6.3) is 
still unknown, although it is almost certain that the list of possibilities is rather small. 
Fortunately, though, the characterization of all complete conformally flat manifolds ad- 
mitting nontrivial solutions to (6.3) was obtained about ten years ago, independently 
by J. Lafontaine [L] and O. Kobayashi [Ko]. The simplification in this case is that 
the Ricci tensor of any conformally flat manifold is always harmonic. Kobayashi and 
Lafontaine proved that beyond the sphere and the product type examples discussed 
above, there are also a collection of examples which are warped products of with 
E^^^ . We shall not write these examples explicitly here, but simply note that, pro- 
vided A; > 2, the solution metrics (S"\A,5f), g e M.Ai are never warped products (for 
topological reasons alone!). We obtain then 

Proposition 6.5. There are no nontrivial solutions to the equation a*{f) for any 
g e M.A, provided k > 2. 

Even when k = 2, only the cylinder M x amongst the Delaunay metrics appears 
on the list of counterexamples. To check this, use the system of equations given in 
Lemma 1.1 of [Ko]. These involve both the warping function for the metric and the 
function /. Elementary manipulations show that the noncylindrical Delaunay metrics 
never satisfy this system. 

Corollary 6.6. Ma is a real analytic Banach manifold in a neighborhood of A4a- 

Proof. It remains only to set up the precise function spaces on which the implicit 
function theorem (see [F] page 239) will be applied, and make a few additional com- 
ments. We use, of course, the weighted Sobolev spaces of symmetric 2-tensors and 
functions, as before. 

The crux of the argument is the fact that when g e A4a, ctg is surjective both on i7| 
and on H^g, and indeed on for any 7 ^ F'! This seems somewhat counterintuitive, 
but follows from the highly over determined nature of a*. In fact, we already know 
that ag has closed range on H^, 7 ^ F' by Proposition 6.4. Its cokernel is identified 
with the kernel of the adjoint, a* which is a map from HZ^ to HZ^~^ ■ By injectivity 
of the symbol of a* we may use elliptic regularity to replace —s here by any positive 
number, e.g. +s. Finally we can invoke the results of Lafontaine and Kobayashi, via 
Proposition 6.5, to conclude that a* has no nuUspace, regardless of the weight —7. 
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Now consider the scalar curvature map 

(6.7) Mg : ff^f (Sym^) x W ^ HUi^j^n^), 

defined in the obvious way. We are including the factor W here so that points of the 
urmoduli space are allowed to have varying neck sizes, but all other perturbations are 
required to decrease exponentially on the ends. To apply the analytic implicit function 
theorem we need to know, first, that the linearization 

(6.8) ag : H'_f{^yra^) ®W ^ H^s 

is surjective, which we have just established, and secondly that (6.7) is a real ana- 
lytic mapping of Hilbert spaces. This latter statement is also straightforward, so this 
completes the proof. 

One unfortunate shortfall of this theorem is that, although Ma provides a thickening 
of M.A, it is not a uniform thickening, i.e. one of fixed 'width' around M.^. However, 
over any compact subset of A^a we can ensure the existence of Ma out to some fixed 
distance by the obvious covering argument. To find this fixed width thickening uni- 
formly on A^A we would need to understand more about the compactification 
If the analogues of the result of Kobayashi and Lafontaine as well as the asymptotics 
Theorem 3.3 and compactness result. Proposition 3.12, were known for non-conformally 
fiat metrics, this would be unnecessary, and many of the results here could be given a 
more satisfactory form. 

Analyticity of Ma and the Generic Slice Theorem 

In addition to its intrinsic nature and relationship to the concerns of this paper, the 
urmoduli space Ma is required for two purposes: to show that M.a itself is analytic, 
and to prove the generic slice theorem. The proofs of these results are very closely 
related, so we develop the preliminaries for them simultaneously. 

Inside the ambient space V = iJf_'^^(Sym^) x W are two analytic (Hilbert) subman- 
ifolds, namely Ma and the conformal class [go]- (By [qq] we always mean the set of all 
metrics in V conformal to the round metric Qq.) Likewise we have all other conformal 
classes in V; a generic such class will be denoted [g']. Notice that all the analytic ma- 
chinery, in particular the Fredholm theory and relative index computations in §4, still 
holds for any metric (corresponding to an element) in V. 

Choose any element g G Ma and let [g'] be the conformal class of g. We shall 
assume that [g'] is very near [go\. The immediate aim is to prove that Ma and [g'] 
intersect almost transversely at 51; in general these submanifolds may not be transverse 
at ^, however they will always form a Fredholm pair there. This means that their 
tangent spaces, E = Tg[g'] and F = T^Ma, which are closed linear subspaces in 
V = TgV = H^_Y © W form a Fredholm pair, i.e. that E D F is finite dimensional, 
and that E + F is closed in V and has finite codimension there. To any Fredholm pair 
one can associate an index, which is the dimension of E (1 F minus the dimension of 
V/{E + F). This number is stable under perturbations of the pair. These matters are 
explained more thoroughly in [K]. 

32 



Lemma 6.9. Ma and [g'] are a Predholm pair at any point g in their intersection. 



Proof. The orthogonal complements of E and F are given by 

E^ = {h:tr,{h) = 0}, 

F^ = {h:h = «*(/), / e Ht_Y}. 

It will suffice to show that EOF and E-^ n F^ are finite dimensional, and that the 
orthogonal projection tt : F — > E-^ has closed range. 

First suppose that he EOF. Then h = f -g with / e H^_Y © W and ag{h) = 0. 
Computing Q;g(/ • g) we find that Lg{f) = 0. Since L is a Predholm operator, there can 
be at most a finite dimensional family of such h in the intersection; the precise dimension 
is the same as the nuUspace of L on On the other hand, if G E-^ fl F-^, then 

h = a*{f) with / e H^'^ and trg(/i) = 0. But trg(Q!*(/)) is, up to a factor, just Lg{f), 
so the dimension of this intersection is the same as the dimension of the nuUspace of L 
on H^^, and is thus finite. 

To conclude the proof, we need to know that the projection tt has closed range (and 
hence is Predholm, by the work above). The map tt is defined, for any h E V, hy 
decomposing h = X ■ g + k, where tr (k) — 0, and setting 7r(/t) = k. Now suppose 
hj is a sequence of elements in F, \\hj\\ = 1 for all j, but kj = 7r{hj) — > in norm. 
Since OLg{hj) = we have (n — l)Lg{\j) = 55kj — {rg,kj). The right hand side of 
this equation certainly goes to zero in norm, so by our Predholm theory for Lg, the 
functions Xj decompose as Xj = + Uj with uj going to zero in norm and in the 
nuUspace. Since Hj -gf is in the nuUspace of ag also, we may subtract it off from hj, and 
get that \\hj\\ 0, contrary to hypothesis. Hence tt restricted to F has closed range, 
and this completes the proof. 

Before stating the next result we need to introduce some notation. Pirst, if [g'] 
is a conformal class on (represented by elements in V, and assumed near to [go]), 
then A1a([5'']) will denote the moduli space of complete metrics on S"^\A with constant 
scalar curvature n{n — 1) in the conformal class [g']. For e > any sufficiently small 
number we also let M.A,e{[g']) denote the subset of MA{[g']) consisting of solution 
metrics g with the Killing norms of all Pohozaev invariants bounded below by e (see 
§7 for a discussion of these norms). This is simply the subset of A^a([^']) consisting 
of solutions with necksizes uniformly bounded away from by e. By Proposition 3.12, 
each A^A,e([(7o]) is compact; then, provided [g'] is sufficiently close to [go], ■MA,e{[9']) 
will also be compact, since it is contained in a compact neighborhood of /Aa^c by virtue 
of the implicit function theorem and the considerations above. 

Generic Slice Theorem 6.11. For any fixed e > 0, the truncated moduli space 
A4A,e{[9']) is a k-dimensional real analytic manifold for all conformal classes [g'] in 
a set of second category and sufficiently near to [gQ\ . 

Proof. The reason for introducing the truncated moduli spaces is, of course, because 
we only know the existence of Ma to a fixed distance away from M.a only over each 
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M.A,e- We shall not comment further on this modification, but simply remark that it 
could be removed provided somewhat more were known about Ma, as discussed at the 
end of the previous subsection. 

This result follows from the Sard-Smale theorem, once we have checked the hy- 
potheses. The basic point, of course, is that M.x{[g']) is the same as Ma n [g']. We 
parametrize the set of conformal classes (not modulo diffeomorphisms!) near to [^o] by 
the linear Hilbert space E-^ = {/c £ H^^^ : tr g{k) = 0} for some fixed g E M.a- Now 
consider the projection map 

(6.12) n : Ma — ^ E^. 

A4a is given by the preimage 11"^ (0). The preimage Il~^{k) in Ma is precisely the 
set M\{[g']) where [g'] = [go + k] is the conformal class corresponding to k. This 
preimage is a smooth, in fact real analytic, manifold provided /c is a regular value of 11. 
By Sard-Smale, once we know that 11 is a Fredholm map (of index k), then the set of 
regular values in E-^ is of second category. Observing that the tangents to the orbits 
of the diffeomorphism group are certainly not of second category, we may conclude 
that the moduli spaces over generic conformal classes near to [^fo] are smooth and k 
dimensional. 

However, the assertion that 11 is Fredholm of index k is contained within the pre- 
ceding Lemma. Indeed, ii g E Ma and F = Tp(Ma), then we proved there that the 
projection tt = II^, : F — E-^ is Fredholm, and its index is the same as the relative 
index of Lg across the weight 0, i.e. equal to k. 

To finish the proof we need to eliminate the possibility that n(MA) is contained 
within a (finite codimensional) submanifold of E-^. Although this would follow from 
knowing that tt is surjective, this is a difficult issue. Instead we appeal to the existence 
theory for the nonlinear equation, generalizing Schoen's basic construction. In fact, 
the construction of solutions given in [S2] may be carried out not only for the standard 
round metric go, but also for generic metrics g' which are small perturbations of go 
compactly supported away from A. The only modification of Schoen's method needed 
to accomplish this is given in [PI], Proposition 1.2. The linearization of this set of 
compactly supported perturbations is clearly dense in i?-*-, hence there is a dense set of 
k G E-^ near the origin for which the preimage Il~^{k) is nonempty. This is sufficient 
to conclude that n(MA) contains a full neighborhood of 0, hence the condition 'regular 
value' for 11 is not the empty one. This completes the proof. 

The reader should note that an alternate possibility for proving smoothness of generic 
slices would be to study the linear operator Lg for any g e Ma near to but not in the 
standard conformal class and show that this operator generically has no kernel. 
The smoothness would then follow by the results of §5, combined with a compactness 
argument. However, proving that Lg has no decaying eigenfunctions would require 
setting up a somewhat elaborate perturbation theory, since this is equivalent to trying 
to perturb point spectrum which is sitting on the end of a band of continuous spectrum. 
Thankfully we have been able to avoid this approach here (although we have had to 

34 



appeal to the rather more diflBcult existence theory for the nonhnear equation instead) . 
It should be noted that the surjectivity of (which we have established) is equivalent 
to this eigenvalue perturbation result. 

The second main result, that M.^ is (locally) a real analytic variety, also follows 
from Lemma 6.9. 

Theorem 6.13. For any g G Mx there exists a hall B in V — i7^"^^(S"'\A; Syrn^) x 
W, a finite dimensional space M , a real analytic variety A in M and a real analytic 
diffeomorphism 



such that $(A1a r\ B) = An B' , where B' is a small hall containing ^{g)- 

The proof reduces to the following abstract result, which is presumably well-known: 

Lemma 6.14. Let V he a Hilhert space, and £ and T two real analytic suhmanifolds, 
the tangent spaces of which at any point of the intersection p & T = S H T form a 
Predholm pair. Then for each point p & T there exists a neighborhood B of p, a finite 
dimensional suhspace M <Z V , a real analytic variety A C M , and a real analytic 
diffeomorphism 

$ : S — >V 

such that $(r r\ B) = Ar\ B' <Z M r\ B' , where B' is a neighborhood of $(p) in M. 

Proof. By initially composing with a real analytic diffeomorphism we may assume that 
a neighborhood p in £ lies in a linear subspace and that p — (In our situation, £ 
is the conformal class and this may be accomplished by a translation.) Let ^ denote 
a 'defining function' for !F in V. By this we mean that ^ is a real analytic map from 
V into another Hilbert space U such that T = ^"-^(0). By the stability of Fredholm 
pairs, any level set ^~^{w) forms a Predholm pair with £, for all sufficiently small 
w & and for all points of intersection near the origin in £. It is evident that this 
statement is equivalent to the assertion that \l/ = : £ U is a Fredholm map. 
The intersection £ H J-' is the same as \E'~^(0). Thus the lemma may be rephrased as 
saying that if is a real analytic Fredholm map between two Hilbert spaces £ and 
[/, then the level set \l/~^(0) is locally equivalent, by an analytic diffeomorphism, to a 
finite dimensional analytic variety. 

A standard implicit function theorem argument proves this assertion. Let A = \I/*|o, 
and set M — ker(A), L — W Q range(A). Then both L and M are finite dimensional. 
Define a new map ^ : £ ® L ^ U hj ^{q, f) ^ *(q) + /. Then 



$ : B 



V 



(0,0) 



{h,cl>) = A{h) + ct>, 



so this differential is obviously surjective. Note that 



ker(** 



(0,0) 



) = {(/i,0):A(/i) + = O}. 
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Since A{h) and (j) lie in orthogonal spaces, this nuUspace is just M = ker(^). The 
analytic implicit function theorem gives the existence of two analytic maps 

k:M^M^c£, £: M ^ L 

and a ball B around (0, 0) C © L such that all zeroes of in this ball lie in the graph 
of the map {k, £) : M ^ © L: 

{{q,f) : ^'(g) + / = 0} nS = {{m + k{m),£{m)) :meMnB}. 

Thus, for m in this ball, '^{m + k{m)) + £{m) vanishes identically. The zeroes of ^ 
in B f] S are then just {m + k{m) : m E B f] M : £{m) = 0}, for these points are 
all zeroes of and are the only such zeroes at which \& also vanishes. If 11 is an 
analytic diffeomorphism of S carrying the graph of k into M, then B fl ^'"^(0) equals 
U~^{MnB) n{m e MnB : £{m) = 0} as desired. This completes the proof of Lemma 
6.14. 

VII. Concluding Remarks and Informed Speculation 

NULLSPACE AND SINGULARITIES OF Ma 

We have not yet discussed whether it is possible to give conditions ensuring that a 
given g e Ma is a 'smooth point,' as defined in §5. It is possible for Ma to be both 
smooth and of the correct dimension near g even if Lg has L^— nuUspace; however, 
absence of this nuUspace is our only criterion for guaranteeing this smoothness. We 
expect that it should be difficult to establish such a criterion in general. It would also 
be quite interesting to understand when Ma is not smooth in the neighborhood of 
some element g. In particular, constructing solutions near which Ma is singular seems 
like another very challenging and important problem. 

On the other hand, the two known constructions for producing points in Ma, those 
of Schoen [S2] and [MPU], yield solutions with explicit geometries. The dipole solutions 
of [MPU] , described below, are manifestly smooth points. These solutions exist only for 
certain configurations A, whereas the ones in [S2] exist for any A. Unfortunately, the 
existence of L^— nuUspace for these latter solutions is less evident, although we expect 
this to be true; outlined below is a strategy to prove this. Note that once the existence 
of one smooth point in any component of Ma is established, the real analyticity of 
the moduli space then shows that almost every element in that component is a smooth 
point. 

We now give a brief description of Schoen's solutions and their construction. These 
solutions are uniformly small perturbations of explicit approximate solutions. Each 
of these approximate solutions is constructed from an 'admissible conformal structure' 
(r, cr); r is an infinite tree with a labeling, a, of strong dilations Gij for each directed 
edge Cij . We assume for simplicity that F has one vertex of order k and all other vertices 
of order 2. The admissibility of the labeling a refers to certain compatibility conditions 
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that the dilations must satisfy. The strength, A^j, of these dilations can be related to 
a measurement of the 'neck sizes' e^j of the approximate solution constructed from 
the data (F, cr). (S"^\A, (/o-) consists of almost spherical regions, corresponding to the 
vertices V of F, joined by small necks, corresponding to the edges, the sizes of which are 
dictated by the e^. There are an infinite number of parameters in the construction of 
these approximate solutions since the admissible conformal structures a can be varied, 
even with F fixed. In particular, one can begin with an initial conformal structure 
o"! , so that the corresponding approximate solution gi consists of one central spherical 
region VLq (the vertex of order k) and k periodic, spherically symmetric ends. But in 
order to find an exact solution, it is necessary to deform the conformal structure, as 
dictated by the Pohozaev balancing condition, and break this symmetry. 

The conformal structure (F, cr) decomposes S"\A into a union of almost spherical 
regions, i eV . Each C is the puUback by a conformal diffeomorphism Fi of 
S"^ of a large region of S". The metric on Vti is constructed so that Fi : (Qj, go) ^ 
(§", ^o) is an isometry off of a small neighborhood of This decomposition gives rise 
to a basic analytic property of the approximate solution metric g^- This is the existence 
of a basis of functions, 77^-, i G and j = l,---,n + l, for an infinite dimensional 
space K (the 'small eigenspace') corresponding to all the spectrum of Lg^ = Ag^ + n 
in a small interval around 0. The functions ry^- have explicit geometric descriptions. 
Each ?7^- has support concentrated on VLi and corresponds to the linear function rjj on 
Fi{Q,i) C Linear functions are of course eigenfunctions for Ag^ with eigenvalue n. 
On the orthogonal complement -ftT-*-, Lg^ is invertible, uniformly in a. As the e^j 0, 
these approximations improve and Lg^ converges, as an operator on L^, to the operator 
Lgg on the disjoint union of spheres indexed by V . 

Writing the solution gf as a perturbation; g — (1 + ry) "-2 g^ and using the conformal 
invariance of the conformal Laplacian, we find that 

Lg=Ag-\-n 

('•1) 4 n+2 

= (1 + ?7)" — Z\g^ +2(1 +?7)-— (V?7, V) +n. 

Now, 7] is small in (by which we mean that Ir^lo can be made arbitrarily small by 
taking all sufficiently small), and estimates from [S2] can be used to show that Vry 
is small in over any fixed compact set, but these facts do not immediately imply the 
existence of an analogous good basis for the 'small eigenspace' K for Lg. Nonetheless, 
we expect this basis to exist, namely that there exist an e > 0, such that if e^j < e for 
all ij then there exists a set of smooth functions, rjij, i E V, j = 1, ■ ■ ■ ,n + 1, which 
satisfy the estimates in Lemma 3.6 of [S2]. The span of these functions is a subspace 
K c L2(§"\A) such that 

||-£'g^7||L2(S"\A) <c(e)||r7||L2(§"\A) ioT r] e K 

||??||l2(S"\a) <c\\Lgr]\\ l2(s^\a) for r] e , 

where c is a constant independent of e, and c(e) tends to as e — > 0. Note that from 
(7.1) this conjecture would be immediate if were finite and small. Unfortunately, 

this is never the case, since the 77 in Schoen's construction is not even in L^. 
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The likeliest method to establish the existence of this space K with its explicit basis 
is to write the solution along each end as a conformal perturbation of the Dclaunay 
metric on the half-cylinder to which it converges at infinity. There is good evidence 
that perturbation determining the solution is uniformly small and even exponentially 
decaying, again provided the neck sizes are sufficiently small. With this information 
in hand, a not very complicated transference procedure would produce the basis for K 
with all necessary estimates. 

If these could be settled, it would then be possible to resolve some basic questions 
about the existence of an L^— nuUspace of L^, for Schoen's solutions. 

Conjecture 7.2. Suppose g e A4a is a solution similar to that constructed in [S2j, 

4 

SO that g = {1 + r])^^-2g^ as above. Assume that A C S"^ does not lie in any round 
hypersphere. Then there exists an e — e(A) > 0, such that if a satisfies eij < e for 
each edge (i.e. all 'neck sizes' are less than e) then Lg has no Lp'—nullspace, and hence 
g is a smooth point of Ma- If ^ does lie in a hyperplane then any element (f) in the 
Lp'—nullspace of Lg is not integrable, i.e. it is not tangent to a path in M.a- 

We provide a brief sketch of our plan for proving this. If is a solution constructed 
from some (r,(T), with all neck sizes sufficiently small, and if e L^, Lgcj) = 0, then 
(j) admits two different decompositions, one into the 'almost linear functions' r]ij con- 
centrated on the almost spherical regions O^, and the other (which is only local along 
each end Eg) into eigenfunctions for the Laplacian along the cross-sectional sphere, 
as in §2. For this latter decomposition we need the (putative) fact that the solution 
stays uniformly close to the Delaunay solution to which it asymptotically converges. 
The former expansion would show that (p is very close to a linear function on each of 
the almost spherical regions. A very important point here is that the linear function 
determined by (p on each of the spherical regions can never vanish. In terms of the 
eigenfunctions of the second decomposition, this first decomposition implies that along 
each end, (p has most of its mass concentrated in the eigenspaces corresponding to the 
zeroeth and first nonzero eigenvalues of the cross-section. Since g is only approximately 
Delaunay, the zeroeth Fourier coefficient (po of 4> is not an exact solution of LqcPq = 0, 
however, the error terms are sufficiently controllable to show that cpo must be quite 
small (again, depending on the size of e). Hence, in fact, most of the mass of is 
contained in the coefficients (/>i, ■ ■ ■ , (/>n, corresponding the the first nonzero eigenvalue 
of the Laplacian on S"^~^ (which has multiplicity n). We call these 'transverse' linear 
functions, since each of the spherical regions where we consider them has a natural axis 
picked out. Thus these arguments show that an solution of Lcp — restricts on 
each spherical region to be almost linear. More specifically, on the innermost sphere 
along each end (the one adjoining the central spherical region), cp actually restricts to 
be approximately transverse linear. But now, on the central sphere, (p is linear and 
near each of the neighborhoods where the various ends are attached, must restrict to a 
transverse linear function relative to the axis determined by that end. This is clearly 
impossible if the points do not lie in any subsphere. 

For the second assertion, concerning the situation when A C we only need to 
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use that 4> is sufficiently close to a linear function on the central sphere. Of necessity, 
this linear function is one which vanishes on the equator determined by A. Hence, if 
were an integrable Jacobi field, we would obtain a family of solutions which were not 
symmetric about this equator, contradicting the reflectional symmetry guaranteed by 
the Alexandrov reflection argument. 

We hope to be able to settle these issues in the near future. There are, of course, a 
number of more refined questions about the singular structure of TWa, beginning with 
the problem, already noted, of constructing solutions near which A^a is singular. 

Coordinates on Ma 

In this section we discuss two ways in which coordinates for M.^ may be given. We 
begin with the one arising from the linear analysis. 

If (7 G AIa is a point for which Lg has no L^— nullspacc, then a neighborhood of g 
in M.\ is parametrized by a small neighborhood of the origin in the bounded nuUspace 
B. Thus, linear coordinates on B yield local coordinates on M.^ near g. The problem 
then is to get precise information about these linear coordinates on B. Unfortunately, 
this seems difficult, in general. We now discuss how one would go about this, and at 
least set the problem up in scattering-theoretic terms. 

By definition, any element (f) in the deficiency subspace W can be expanded on each 
end as a combination of the model solutions 0i, 02 plus an exponentially decreasing 
error. For clarity, we label these model solutions on Ej with a superscript [j). Thus 

(7.3) W3(l)^ aj(f)[^^ + hj(j)i^ on Ej. 

So any (f) & W determines a map 

S:W — ^ R^'^ 
01 — )• (ai,6i,...,afe,6fe) 

By the definition of W, S is an isomorphism. Now suppose is in the bounded nuUspace 
B. 

Lemma 7.5. Under the hypothesis (5.2), the linear map 

C = S\^:B — > 

is infective. Its image Sa is a k-dimensional subspace which is Lagrangian with respect 
to the natural symplectic form Yl^=i ^'^j ^ ^^j ■ 

Proof. C is obviously linear, and it is injective, because otherwise there would exist an 
element (p E B with all its asymptotic coefficients aj, bj vanishing. By the results of the 
last section, such a would lie in L^, which contradicts the hypothesis (5.2). 

That ^A is Lagrangian is a consequence of Green's theorem. There is an analogous 
result for manifolds with asymptotically cylindrical ends proved in [Me]. Since we 
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(7.4) 



already know that S\ is /c-dimensional, it suffices to prove that it is isotropic. For this, 

let tp, and set C{ij) = (ai, 61, . . . , a^, 6^), C{x} = ("i, A, • • • , Pk)- Using the 

coordinate tj along the end Ej, let Ma denote the region in S"\A where each tj < A. 
Then, since Lip — Lx — 0, 

amx - ^{Lx)) dV, = / (^x - da„ 
Ma JdMA 

where u is the unit vector-field normal to OMa- Substitute the expansions for ip and 
X at each end, and drop all terms which decrease exponentially as ^ — > 00 in this 
expression. This means that we can replace the metric g on Ej by the Delaunay metric 

g^- there, and also replace v by u^f^^"^ '^^^tj and dag by ui^^ Taking the limit 

as ^4 — > 00 we get 

(7.6) = tiaA-a,Ml_^ rf,,. 

Here da is the standard volume measure on If we set Aj equal to the integrand for 

the j^^ end, then it suffices to show that these constants Aj are nonzero and independent 

of j. However, (p^^^ and (p2 '^ are independent solutions of the ODE Lgcf) = 0, and the 
expression in parentheses in the integrand in (7.6) is just the Wronskian for this ODE. 
This Wronskian may be written explicitly as some nonzero multiple of 

exp I -2 / ^ 1 = u;' 

This shows that the integrand, hence the integral itself, is a nonvanishing constant, 
independent of j. 

The relevance of 5'a to the description of coordinates on Ma near g is explained as 
follows. Suppose first, for simplicity, that there exists an element (/> e i3 for which all 
ttj, bj — 0, j > 1, but ai = 1, 61 = 0. If g{s) is a curve in moduli space tangent to this 
0, then the Delaunay parameters ej on each of the ends would remain fixed, at least 
infinitesimally, along this curve. In fact, the metrics g{s) (or the underlying curve of 
solutions Us of the PDE) would be strongly asymptotic to the initial metric g on Ej, 
j > 1. On El however, g{s) would be strongly asymptotic to a translation of g, or 
alternately, of its model Delaunay metric there. If, on the other hand, ai = 0, bi = 1 
and all other aj, bj = 0, then again each g{s) would be strongly asymptotic to g on all 
ends except Ei; on that end the Delaimay parameter for the model metric would be 
changing, but the model Delaunay metric would not be translating. Of course, neither 
of these types of elements of B need exist. 

If is a general element of B, and g{s) is a curve in A4a through g tangent to (p, then 
(p describes infinitesimal changes in translation and Delaunay parameters along each 
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end. The most precise description of coordinates on A^a would relate the proportions 
of these various changes on each end to one another for all directions (p G TgM.\. This 
is equivalent to describing the coefficients aj,hj for each (p & B, and this, in turn, is 
equivalent to describing the Lagrangian subspace Sa. An explicit description of Sa 
requires a better understanding of the Pohozaev balancing condition (3.9). 

For the sake of illustration, let us examine all this for A42, the moduli space when 
k = 2. Of course, we have already provided a complete description of this space, but it 
provides a concrete example of this description. We let the background metric be any 
Delaunay metric ge with e < u. Now is 2 dimensional, consisting of the elements (pi 
and 02- Hence 

5A = {al = a2,6l = 62}cR^ 

which is Lagrangian, as expected. The two natural curves g{s) emanating from any 
g e M.2 are simply the families $i and $2 defined in (2.14). 

It may be possible to give a geometric description of parameters on Ma by consid- 
ering the Killing norms of the Pohozaev invariants. Recall that these invariants are 
elements of o*(r2, + l, 1), the dual of the Lie algebra of the conformal group. Any element 
X e o(n + 1, 1) may be uniquely decomposed as X = Xq + w, where Xq e o(n + 1) 
and w e M"+^ The Killing form 

B :o{n + 1, 1) X o(n + 1, 1) ^ R 

is the nondegenerate symmetric quadratic form given by 

B{X, X) = ^Tr{XoXo) + w • w. 

Thus, B{-, •) is positive definite on and negative definite on o(n-l-l) (with respect 

to this decomposition). Moreover, the adjoint representation preserves B in the sense 
that B{Ad{F){X),Ad{F){X)) = B{X,X), for aU F e 0(n+ 1,1) and X,X e o{n + 
1,1). 

Since B is nondegenerate it provides an identification between o(n + 1, 1) and its 
dual space. Thus we may use the Killing form to identify, for the Pohozaev invariants 
Vi, . . . ,Vk of a, metric g G A^a, corresponding elements V^^ . ■ ■ ,V'^ G o(n + 1, 1). The 
squared killing norms of the Pohozaev invariant Vi, is then defined to be 

\\r^\\^^B{viv'^. 

These k numbers are conformal invariants. By this we mean the following. If F G 
0{n + 1, 1) is a conformal diffeomorphism of and g G Ma then F'*'{g) G A^^-i(a)- 
The Pohozaev invariants of g and F*{g) do not coincide but transform via the co-adjoint 
representation (see [S2] and [KKMS]). Since B is invariant under this representation 
this implies that the squared killing norms of the Pohozaev invariants of g and F*{g) 
are the same. 
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Schoen has suggested the foUowing way to obtain parameters on Ma- Let A = 
{pi, . . . be a balanced singular set. One can then try to produce, for some e > 0, a 
1-parameter family of solutions, gt G M.x, t G (0, e), such that the asymptotic neck sizes 
ei, . . . , e/c for gt are all equal to t. To realize the other {k — l)-parameters, choose any 
A;-tuple of numbers, (ai, . . . , Ofc), close to 1 and normalized so that ai = 1. Then there 
is a conformal diffeomorphism F taking A to A' = {p[, . . . so that Yl'j=i ^iP'j — 0- 
As before, there is now a 1-parameter family of solutions g'^, on \ A' with asymptotic 
neck sizes given by (ait, . . . ,akt). Thus F*{g[) G A^a is another 1-parameter family 
of solutions. This exhibits the fc-parameters as (t, a2, . . . , a^). It should be possible to 
phrase this in terms of the squared Killing norms of the Pohozaev invariants described 
above. 

DiPOLE Solutions 

We briefly describe a new construction of solutions [MPU] for certain special config- 
urations A. 

Because the Delaunay solutions on the cylinder are nondegenerate, in the sense that 
the linearized operators have no L^— nuUspace, it is relatively easy to graft them 
onto known nondegenerate solutions. In this way, it is possible to construct an explicit 
space of solutions when A has 2k points which are grouped in pairs cr^, i = 1, . . ./c, 
where the distance between points in each pair is very much less than the distance 
between pairs. We state the result when k = 2 and the two pairs are assumed to be 
located at opposite sides of 

Theorem 7.7 ([MPU]). Given e > there exists an rj > 0, such that for any A = 
{pi,P2,P3,P'i}, with d{pi,p2) < rj, d{p3,p4) < rj, and d{pi,p3) > 1, there exists a four 
parameter non- degenerate moduli space of solutions to the singular Yamabe problem 
with neck-sizes at the points pi all greater than e. Moreover, the map from to M'^ 
sending a solution to the translation and neck parameters of one point in each of the 
pairs, say pi and ps, is an isomorphism onto an open set. 

These solutions can be described by making a conformal change of the Delaunay 
solutions to R"\{gi, 52} and R"'\{q'3, q'4}. Simply adding the two solutions gives a result 
which is a very close approximation to an exact solution which may be constructed 
using the implicit function theorem. It is necessary in this construction to make a 
small change in the neck size and translation parameters of one point in each pair to 
make the linear operators here surjective on the relevant spaces. Unfortunately, no 
uniqueness theorem has been proved, even in this restricted setting. 

It is of interest to note that since the parameters on the other point of each pair 
may be held fixed, this construction produces solutions with asymptotically cylindrical 
ends. 

The Structure of the Boundary OMa 

There are a number of interesting questions about the structure of the compactifi- 
cation M.a- Recall from §3 that this compactification is obtained by adding certain 
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solutions g e M.\i to AIa, for certain subsets A' C A. There are two basic problems. 
The first is to determine which subsets A' have elements occurring in dM.;^. This is 
again intimately connected to the balancing condition (3.9). There are few examples 
where we can say anything explicit about this. The case k — 2 was already discussed in 
§3. The existence of dipole solutions [MPU] may give additional information concerning 
this problem when k is even. 

The other fundamental problem is to determine how much of M.a' is contained in 
dM-x when at least one point g e M.k' lies in dM.\. More specifically, we propose the 
following: 

Conjecture 7.8. Suppose that g G M.a' H OMa- Then the entire component of M.a' 
containing g also lies in dM.A- 

It is easy to show, using the compactness theorem of [P2], that the set of points g' , 
in the same component of M.a' as g and which also lie in dM-A-, is closed. One method 
to prove that this set is open would be to generalize the construction of dipole solutions 
to a more general grafting procedure. This would prove the conjecture. 
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